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ABSTRACT 



When the neutral interstehar medium is exposed to EUV and soft X-ray radiation, the argon 
atoms in it are far more susceptible to being ionized than the hydrogen atoms. We make use of 
this fact to determine the level of ionization in the nearby, warm, neutral medium (WNM). By 
analyzing FUSE observations of ultraviolet spectra of 44 hot subdwarf stars a few hundred pc 
away from the Sun, we can compare column densities of Ar I to those of O I, where the relative 
ionization of oxygen can be used as a proxy for that of hydrogen. The measured deficiency 
[Ar I/O I] = —0.427 ±0.11 dex below the expectation for a fully neutral medium implies that the 
electron density n{e) w 0.04 cm~'^ if n(H) = 0.5 cm~'^. This amount of ionization is considerably 
larger than what we expect from primary photoionizations resulting from cosmic rays, the diffuse 
X-ray background, and X-ray emitting sources within the medium, along with the additional 
ionizations caused by energetic secondary photoelectrons. Auger electrons, and photons from 
helium recombinations. We favor an explanation that bursts of radiation created by previous, 
nearby supernova remnants that have faded by now may have elevated the ionization, and the 
gas has not yet recombined to a quiescent level. A different alternative is that the low energy 
portion of the soft X-ray background is poorly shielded by the H I because it is frothy and has 
internal pockets of very hot. X-ray emitting gases. 

Subject headings: ISM: atoms - ISM: general - ISM: lines and bands - local interstellar matter - tech- 
niques: spectroscopic - ultraviolet: ISM 



In large part, we are aware of the principal pro- 
cesses that create free electrons in otherwise com- 
pletely neutral parts of the interstellar medium 
(ISM) in the disk of our Galaxy. However, for 
some of the contributing factors, ones that either 
enhance or diminish the relative ionization of the 
gas, there is a need to validate our understand- 
ing of their strengths. Most of these processes 
are well understood qualitatively, but one major 
quantitative uncertainty is the effectiveness of ex- 
treme ultraviolet (EUV) and soft X-ray photons 
in ionizing the gas. Two factors contribute to this 
uncertainty: one is the difficulty in measuring the 
fluxes of diffuse photons with energies above the 
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ionization potential of hydrogen (13.6 eV) but be- 
low about 100 eV, and the other is an overall as- 
sessment of how well these photons can penetrate 
the neutral regions, which depends on the poros- 
ity of the gas structures and the distribution of 
radiation sources. 



Our ultimate goal is not only to understand 
these processes better, but also to obtain an es- 
timate for the fractional amount of the gas that is 
in an ionized state. This information is relevant 
to gauging the strength of heating of the gas due 
to the photoelectric effect from dust grains irra- 
diated by starlight, since the grain charge, which 
regulates its rate, depends on the electron density 
(Weingartner & Draine 2001b). While generally 
considered to be less important than the photo- 



electric effect, other means of creating thermal en- 
ergy, such as the heating caused by secondary elec- 
trons from cosmic ray and X-ray ionizations and 
the dissipation of Alfven waves and magnetosonic 
turbulence (Kulsrud & Pearce 1969; Mclvor 1977; 
Spangler 1991; Mintcr & Spangler 1997; Lerche et 
al. 2007), depend on the relative fractions of elec- 
trons. Cooling of the gas through the collisional 
excitation of the fine-structure levels of C"^, Si"*" 
and Fe+ or, at temperatures approaching 10^ K 
various metastable levels and the La transition of 
hydrogen, are likewise governed by the degree of 
partial ionization (Dalgarno & McCray 1972). In 
the following paragraphs, we consider three differ- 
ent pathways for creating small amounts of ioniza- 
tion in the mostly neutral ISM. 

Except for the interiors of dense clouds where 
there is significant extinction, all of the space 
in the disk of our Galaxy is exposed to ultra- 
violet starlight photons that are capable of ion- 
izing atoms that have first ionization potentials 
less than that of hydrogen (13.6 eV). Since the 
recombination rates of the ions are slow relative 
to their ionization rates, the concentrations of the 
ionized states of these atoms are strongly domi- 
nant. Thus, it is a simple matter to add together 
the contributions from various elements that are 
able to supply free electrons. The only uncertainty 
here is an accounting of the strengths of depletions 
from the gas phase caused by these elements con- 
densing into solid form onto dust grains. These 
strengths vary collectively for all of the elements 
from one region to the next. If we take such vari- 
ations into account (Jenkins 2009), we can state 
that most of the gas will have free electron con- 
tributions that should be somewhere in the range 
n(M+) = 0.8 - 1.7 X 10-'*n(Htot), where n{M+) 
is the number density of heavy elements that are 
capable of being ionizec^ and n(IItot) is the total 
density of hydrogen in both neutral and ionized 
forms. 

As cosmic ray particles collide with gas atoms 
in the Galaxy, they heat and ionize the ISM. We 



^This range was computed from the solar abundances and 
representative values of the gas fractions [Xgas/H]i and 
[Xgas/H]o listed in Jenkins (2009) for elements that can 
be ionized by starlight, except that the gas-phase abun- 
dance of C was lowered by a factor of 0.43, in accord with 
the recommendation by Sofia et al. (2011) that earlier de- 
terminations of A'^(C II) were systematically too high by a 
factor of 1/0.43. 



are unable to observe the flux of the lowest en- 
ergy particles because we are shielded from them 
by the heliospheric magnetic field, and extrapola- 
tions from the observed higher energy flux distri- 
butions are uncertain (Spitzer & Tomasko 1968). 
Nevertheless, from measurements of the relative 
abundances of trace molecular species, the cosmic 
ray ionization rates ^CR ~ 0.5 — 9x 10~^^s~^ seems 
to be the most plausible range for the general ISM 
(Wagenblast & WiUiams 1996; Liszt 2003; Indriolo 
et al. 2007; Neufeld et al. 2010; Indriolo & Mc- 
Call 2012), although details in the chemical mod- 
els may introduce some uncertainty [cf. Le Petit 
et al. (2004) and Shaw et al. (2008)]. The chem- 
ical models of Bayet et al. (2011) suggest that 
in some particularly active regions the ionization 
rates may increase to CcR > 1 x 10~^^s^^. 

We now consider a third mechanism for ioniz- 
ing the gas, one that is harder to quantify than 
the other two. EUV and soft X-ray radiation 
can ionize atoms, through both the action of pri- 
mary photoionizations and by creating a cascade 
of energetic, secondary photoelectrons that can 
coUisionally ionize other atoms. Estimates of the 
effectiveness of these agents are difficult to syn- 
thesize, since there are many complicating fac- 
tors. At high energies, much of the radiation 
arises from cooling, very hot {T > 10® K) gas 
coming from recently shocked regions within the 
disk and halo of our Galaxy. For energies slightly 
above about 100 eV, the photons can survive a 
journey through the neutral medium up to about 
Af(H°) = 2 X 10^^ cm^^, and this penetration 
depth progressively increases with energy. Sup- 
plementing this ionizing radiation is that coming 
from several kinds of sources that are embedded 
within the neutral medium. These include stars 
over a wide range of spectral types on the main 
sequence, active X-ray binaries, and the plentiful, 
but faint white dwarf (WD) stars. 

2. STRATEGY OF THE INVESTIGA- 
TION 

The objective of the current study is to use 
specialized observations to help resolve the un- 
certainties that were mentioned above and gain 
a quantitative insight on how effectively the neu- 
tral regions are partially photoionized by EUV 
and X-ray radiation. We do this by repeating a 
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method developed by Sofia & Jenkins (1998, here- 
after SJ98), who examined interstellar UV absorp- 
tion features that could be seen in the spectra of 
background stars so that they could compare the 
abundance of the neutral form of argon, which is 
highly susceptible to being photoionized, to that of 
hydrogen, which has an ionization potential very 
close to that of argon but with a markedly lower 
ionization cross section. We propose that it is 
safe to assume that the abundance of argon (both 
neutral and ionized) relative to that of hydrogen 
should be equal to the solar ratio. For instance, 
SJ98 presented arguments that support the prin- 
ciple that argon in the gas phase is not likely to 
be depleted by being incorporated into an atomic 
matrix within interstellar dust grains. We will re- 
inforce this idea with some indirect observational 
evidence in Section 14.2.11 Thus, we operate on 
the principle that any deficiency of neutral argon 
(Ar I) below our expectation for the amount of H I 
that is present may be considered to arise from the 
conversion of Ar to its ionized form, which is in- 
visible. 

The investigation by S J98 covered an extremely 
limited number of target stars observed with the 
Interstellar Medium Absorption Profile Spectro- 
graph (IMAPS) (Jenkins et al. 1996). Later, 
Jenkins et al. (2000b) and Lehner et al. (2003) 
reported on observations of absorption features of 
Ar I observed by the Far Ultraviolet Spectroscopic 
Explorer (FUSE) (Moos et al. 2000; Sahnow et al. 
2000) toward collections of WD stars inside and 
slightly beyond the edge of the Local BubbleU in 
order to infer ionization conditions of clouds sub- 
jected to the characteristic radiation field in our 
immediate neighborhood. Our current goal is to 
extend our reach well beyond the stars surveyed 
in these two studies, again by using FUSE spec- 
tra, so that we can sample regions of more typical 
densities well outside the Local Bubble. We do 
this by downloading from the Mikulski Archive for 
Space Telescopes (MAST) at the Space Telescope 
Science Institute a large collection of FUSE spec- 
tra of hot subdwarf stars that are situated several 



^The Local Bubble is an irregularly shaped region with an 
unusually low average density with a radius of about 80 pc 
that is approximately centered on the Sun (Vergely et al. 
2010; Welsh et al. 2010; Reis et al. 2011). It contains small, 
partly ionized clouds immersed in a much lower density 
medium (Redfield 2006; Redfield & Linsky 2008). 



hundred pc away from us, well beyond the bound- 
ary of the Local Bubble. 

Our ultimate objective is to compare the neu- 
tral fractions of Ar and H, as had been done in 
past investigations. However, a conventional ap- 
proach of simply deriving the two column densi- 
ties and dividing one by the other is not easily 
achievable with the data in this survey for two 
reasons. First, it is difficult to measure iV(Ar I) 
because the lines are saturated, but only moder- 
ately so, and recorded at low resolution. Second, 
the amount of H I on a sight line can usually be de- 
termined from the damping wings of La, but when 
the interstellar column density is low and there is 
a significant stellar La absorption, one must know 
the star's effective temperature and surface grav- 
ity and then create a model for the stellar feature, 
against which the interstellar feature is superim- 
posed. Also, as emphasized by Sofia et al. (2011) 
measurements of A^(H I) using the damping wings 
of La can give a misleading outcome if one does 
not know about and correct for the effects of the 
velocity structures of the gas. This problem is 
probably most severe for low column density cases 
in the present study. 

We can overcome the difficulties mentioned 
above if we replace H with O as the comparison el- 
ement. In the wavelength coverage of FUSE there 
are a large number of O I absorption features, 
and these lines cover a broad range of transition 
probabilities. Most important, the strengths of 
the O I features are comparable to the one avail- 
able feature of Ar lH By a simple comparison of 
the strength of this Ar I line to those of O I, we 
can determine 7V(Ar I)/A^(0 I). We describe this 
process in more detail in Section 221 

The ionization fraction of O is strongly locked 
to that of H through a strong charge exchange 
reactiorlfl (Field & Steigman 1971; Chambaud et 



•^There are two transitions of Ar I in the FUSE wavelength 
coverage. We can use only the one at 1048.220 A because 
the 1066.660 A line has interference from a pair of strong 
stellar Si IV lines at nearly the same wavelength. 

■^We add a caution that deviations from a nearly one-to-one 
relationship for the ionization fractions of O and H can oc- 
cur at low temperatures because the ionization potential of 
O is slightly higher than that of H (AE/k = 229 K) for the 
lowest fine-structure level in the O I ground state. How- 
ever, since the ionization fractions are small and most of the 
gas we are considering is at temperatures much higher than 
AE/k, this deviation is generally small enough to ignore. 
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al. 1980; Standi et al. 1999). Thus we can use 
O as a proxy for H. The only shortcoming of this 
tactic is that O can be depleted in the ISM, but the 
depletion factors are not very large in the regimes 
of low densities considered here (Cart ledge et al. 
2004, 2008; Jenkins 2009). 

Section [3] of this paper describes the selection 
of archival FUSE spectra and how they were pro- 
cessed to yield useful presentations of the ab- 
sorption features for measurements of equivalent 
widths. Our method of interpreting the spectra 
to yield the ionization of Ar relative to that of 
O is presented in Section SI Section [5] contains a 
short digression on how we verify that the target 
stars are beyond the edge of the Local Bubble. In 
Section [6] we outline the basic equations that take 
into account the processes that influence the par- 
tial ionizations of Ar, relative to those of H and O. 
The equations presented in this section are virtu- 
ally identical to those outlined by SJ98, but with 
some new refinements (i.e., a few reactions that 
were not included earlier). We consider the cre- 
ation of free electrons from the starlight ioniza- 
tion of heavy elements and the effects from cosmic 
rays as processes that are mostly understood and 
already accounted for, and view the actions arising 
from EUV and X-ray ionizations as the principal 
unknowns whose strengths are to be determined. 
In Section [7] we make a prediction for the degree 
of ionization produced by known sources of radia- 
tion, but find that in order to satisfy the general 
outcomes for our measured ratios of Ar I to O I, 
an extraordinarily low volume density of hydro- 
gen rt(Htot) is required. In order to obtain the 
same results for higher densities, we must propose 
a means of achieving higher levels of ionization. In 
Section[8]we propose two possibilities: (1) there is 
a large residual ionization left over from effects 
of radiation emitted by nearby, but now extinct 
supernova remnants (SNRs) over the past several 
Myr or (2) the neutral medium is porous enough 
to allow external, low-energy photons to penetrate 
the gas with less than the expected amount of at- 
tenuation. Section |9] presents an overview of the 
implications of our results on an assortment of 
physical processes and various other kinds of ob- 
servations that depend on electron fractions in the 



Nevertheless, we performed explicit calculations of the O 
and H ionization fractions, as described later in Section [6.3l 



diffuse, neutral medium. The paper ends with a 
summary of the main conclusions (Section [10]). 

Appendices to this paper give descriptions of 
various atomic processes that were incorporated 
into the calculations, but at a level of detail that 
most readers may wish to ignore. A general sec- 
tion on the ionizations arising from secondary elec- 
trons (Appendix |A| is broken into two subsections: 
one treats the effects from electrons liberated by 
the ionizations of H and He (Section lA.ip . while 
the second one discusses primary and Auger elec- 
trons created by the inner shell X-ray ionizations 
of heavy elements (Section lA.2p . AppendixlBl gives 
the equations for evaluating the effects from a mul- 
titude of different kinds of ionizing photons that 
arise from the recombination of IIe++ and He+ 
ions with free electrons. Finally, Appendix [C] de- 
scribes how we can estimate the rates of cosmic 
ray ionization of H°, He" and Ar° from the ob- 
served rates that apply to molecular hydrogen in 
dense clouds. 

3. OBSERVATIONS 
3.1. Target Selection 

Our objective was to make use of target 
stars that represented intermediate cases between 
nearby WD stars, whose sight lines are entirely 
or heavily influenced by conditions in the Local 
Bubble, and the much more distant hot main- 
sequence, giant or supergiant stars that can create 
their own enhanced ionizations in atypical con- 
centrations of gas associated with their formation. 
Hot subdwarf stars represent a class of objects 
that fall into this intermediate category. They 
have distances that are of order a few hundred 
parsecs from us, which reduces the contribution 
from material in the Local Bubble to a very minor 
level. Since they are old, they have had adequate 
time to escape from their progenitorial gas clouds, 
and thus their locations are essentially random 
and should show no preference for dense gas com- 
plexes. They have the additional advantage that 
they are bright enough to yield good quality spec- 
tra, but they are not so bright that they exceed 
the maximum allowed count-rate levels for FUSE. 

In an initial screening of prospective targets, 
we examined the quick-look plots of all stars clas- 
sified as sdO and sdB spectral types in the MAST 
archive of FUSE data. In this step, we rejected 
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all spectra that either seemed to show very strong 
stellar spectral features (8 stars) or that had an 
inadequate signal-to-noise ratio (S/N) at wave- 
lengths in the vicinity of 920 A (70 stars), which 
is where most of the O I lines are situated. A 
few further rejections were made after the spec- 
tra were downloaded and found to have observ- 
ing anomalies (an extraordinarily large number of 
missing observations caused by channel misalign- 
ments: 3 stars), strong stellar features that were 
not evident in the quick-look plots (3 stars), or 
molecular hydrogen lines that were strong enough 
to seriously compromise the lines that we wanted 
to measure (2 stars). The lack of stars with ex- 
ceptionally strong H2 features helped to eliminate 
sight lines that penetrate dense, cool gas clouds. 

3.2. Creation of the Spectra 

All of the downloads of the calibrated FUSE 

data from MAST occurred well after the final 
pipeline reductions were performed for the archive 
with CalFUSE version 3.2.3 (Dixon et al. 2007). 
For every target, wc accepted data from all of the 
available observing sessions (identified by unique 
archive root names) but rejected any subexposures 
that had an extraordinarily low count level caused 
by a channel misalignment during the observation. 
We used exposures obtained during both orbital 
day and night. Normally, one must be cautious 
about observations of features for either O I or 
N I because they can be filled in by diffuse telluric 
emission lines during daytime observing. However 
the O I transition strengths considered here are so 
weak that the telluric contributions are insignifi- 
cant. 

All subexposures and spectral channels that 
passed our initial screening were coadded with 
weight factors based on the inverse squares of their 
respective values of S/N for intensities smoothed 
over a wavelength interval of 0.12 A (or 9 in- 
dependent spectral elements this ensiircs that 
weights are not strongly influenced by random 
noise excursions). Before this coaddition took 
place, wc examined some strong interstellar fea- 
tures and aligned the individual spectra in wave- 
length against a preliminary coaddition with no 
wavelength shifts. This process enabled us to 
virtually eliminate any degradation in resolution 
caused by drifts of the spectra in the wavelength 
direction from one subexposure to the next. How- 



ever, there can still be overall small systematic er- 
rors in radial velocity of about 10 km s~^; see Ap- 
pendix A of Bowen et al. (2008). In a few cases, 
the spectral S/N values were too low to allow such 
shifts to be made with much confidence, even after 
the intensities were smoothed with a median filter 
for viewing. Such spectra were combined without 
any shifts. For every target, two combined spectra 
were created: one was made up with shifts appro- 
priate to the spectral region covering the Ar I line 
at 1048 A, while the other had differentials that 
were optimized for the wavelengths that covered 
the weakest O I lines near 920 A. 

4. ANALYSIS 

4.1. Equivalent Widths and Their Errors 

We measured equivalent widths of the Ar I 
and O I lines by integrating intensity deficits be- 
low best-fitting Legendre polynomials for the con- 
tinua defined from intensities at locations some- 
what removed from the features. Special precau- 
tions were made to account for various sources of 
error, which are important for later analysis stages 
that assign relative weights to different measure- 
ments and also for the estimates of the ultimate 
errors in the results. First, we accounted for the 
direct effect that random count-rate variations can 
have on the equivalent width outcome (Jenkins et 
al. 1973). Next, the weakest lines are subject to 
uncertainties arising from improper definitions of 
the continua. To construct the probable errors, we 
evaluated the expected formal errors in the poly- 
nomial coefficients, as described by Sembach & 
Savage (1992), and then we multiplied them by 
2 in order to make an approximate allowance for 
additional uncertainties caused by the arbitrari- 
ness in selecting the most appropriate polynomial 
order. To find the effects of these continuum er- 
rors on our measurements, the equivalent widths 
were re-evaluated using the probable excursions of 
the continua on either side of the preferred ones. 
Errors in the background subtraction in the FUSE 
data processing are small compared to the other 
errors. 
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Table 1 
Equivalent Widths'* 



Star Name 



Ar I" 
1048.220 
2.440 



O 1''''= 



(1) (2) 

2MASS15265306 

+7941307.. 77 ±18 

AA Dor 60 ± 10 

AGK+81 266. . . 76 ± 19 

BD+18 2647... 43 ± 18 

BD+25 4655 . . . 33 ± 18 

BD+28 4211 ... 31 ± 8 

BD+37 442 166 ± 7 

BD+39 3226 ... 54 ± 3 

CPD-31 1701. . 16 ± 8 

CPD-71D172.. 34 ±10 

EC11481-2303. 132 ± 16 

Feige 34 36 ± 17 

HD 113001 115 ± 20 

JL 119 112 ± 5 

JL 25 77 ± 19 

JL 9 100 ±10 

LB 1566 57 ±6 

LB 1766 67 ±16 

LB 3241 56 ± 9 

LS 1275 62 ± 15 

LSE 234 115 ±12 

LSE 259 149 ±6 

LSE 263 119 ±17 

LSE 44 93 ±12 

LSII+18 9 91 ±15 

LSII+22 21 45 ±12 

LSIV+10 9 179 ± 12 

LSS 1362 88 ± 3 

MCT 2005 

-5112 103 ±8 

MCT 2048 

-4504'' 122 ± 12 

NGC6905 

star" 147 ± 8 

PG0919+272 ... 55 ± 27 

PG0952+519... 52 ±3 

PG1032+406 . . . -21 ± 27 

PG1051+501 . . . 165 ± 15 

PG1230+068... 102 ±15 

PG1544+488 ... 84 ± 3 



919.917 
-0.788 
11.6 
(3) 



922.200 

-0.645 
15.3 
(4) 



925.446 
-0.484 
20.6 
(5) 



916.S15 
-0.362 

25.0 

(6) 



930.257 


919.658 


921.857 


924.950 


950.885 


976.448 


94S.686 


971.738 


-0.301 


-0.137 


-0.001 


0.155 


0.176 


0.509 


0.778 


1.052 


27.8 


34.5 


40.4 


46.9 


49.0 








(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


24 ± 22 


41 ± 22 


39 ± 19 


84 ± 11 


59 ± 15 


97 ± 12 


111 ± 11 


(120 ± 12) 


75 ± 31 


51 ± 37 




93 ± 13 


73 ± 21 


(121 ± 19) 




(225 ± 9) 


60 ± 19 


76 ± 15 


87 ± 18 


99 ± 14 


122 ± 14 


(116 ± 18) 


(102 ± 15) 


(125 ± 14) 


26 ± 16 


39 ± 18 




78 ± 11 


68 ± 11 


(63 ± 13) 


(131 ± 12) 


(159 ± 11) 


22 ± 11 


21 ± 12 




53 ± 11 


41 ± 12 


64 ± 13 


77 ± 12 


(75 ± 11) 


37 ± 4 


34 ± 4 


35 ± 4 


42 ± 3 


51 ± 3 


(42 ± 4) 


(83 ± 4) 


(62 ± 3) 


95 ± 24 




166 ± 17 


128 ± 16 










48 ± 10 


47 ± 8 


38 ± 8 


52 ± 7 


59 ± 7 


54 ± 7 


69 ± 8 


(73 ± 7) 


20 ± 11 


39 ± 6 




45 ± 4 


(28 ± 5) 


(45 ± 7) 


(68 ± 5) 


(91 ± 4) 


57 ± 13 


47 ± 15 


36 ± 18 


67 ± 12 


67 ± 13 


(61 ± 13) 


(93 ± 13) 


(83 ± 14) 


92 ± 16 


99 ± 17 


80 ± 34 


145 ± 15 


121 ± 17 


153 ± 17 


146 ± 16 


177 ± 16 


28 ± 14 


38 ± 13 




42 ± 13 


46 ± 14 


41 ± 21 


86 ± 19 




(53 ± 15) 


63 ± 10 




96 ± 10 


93 ± 12 


100 ± 13 


125 ± 11 


128 ± 10 


103 ± 15 


83 ± 25 


80 ± 23 


138 ± 7 


111 ± 19 


(169 ± 7) 


(180 ± 6) 


(174 ± 7) 


72 ± 22 


29 ± 17 




69 ± 17 


94 ± 17 




111 ± 18 


106 ± 17 


95 ± 16 


87 ± 14 




106 ± 9 


116 ±9 




129 ± 11 


(142 ± 12) 


56 ± 14 


48 ± 15 




81 ± 11 


73 ± 19 


(98 ± 11) 


(108 ± 14) 


(111 ± 11) 


58± 14 


38 ± 14 




98 ± 11 




120 ± 11 


(148 ± 10) 


(168 ± 10) 


53 ± 8 


52 ± 8 


97 ± 8 


82 ± 7 


106 ± 7 


(112 ± 7) 


(123 ± 7) 


(128 ± 7) 


65 ± 13 


46 ± 10 


(99 ± 13) 


.58 ± 9 


74 ±10 


80 ± 10 


89 ± 11 


(89 ± 9) 


94 ± 16 


85 ± 18 


123 ± 17 


107 ± 1.5 


122 ± 1.5 




154 ± 18 


1.30 ± 15 


134 ± 20 






123 ±13 


156 ± 13 


■ ■ ■ 






80 ± 15 


76 ± 14 


■ ■ ■ 


91 ± 13 


76 ± 13 


95 ± 17 


139 ± 14 


159 ± 12 


78 ± 15 


61 ± 15 


94 ± 12 


92 ± 9 


132 ± 10 




(^140 ± iUJ 


±1-1) 


73 ± 23 


72 ± 21 




78 ± 21 


125 ± 21 




(119 ± 21J 


52 ± 13 


44 ± 11 


(116 ± 8) 


63 ± 8 


74 ± 9 


(91 ± 8) 


(100 ± 9) 


(119 ± 8) 


106 ± 18 


88 ± 12 




152 ± 9 


170 ± 10 


231 ± 9 




(192 ± 10) 


78 ± 18 


64 ±20 


(65 ± 8) 


(68 ± 10) 


134 ± 7 








75 ± 18 






113 ± 9 


135 ± 12 


(168 ± 15) 


(150 ± 13) 


(124 ± 18) 


61 ± 59 




90 ± 28 


118 ± 14 


117 ± 17 


149 ± 14 


174 ± 16 


(172 ± 15) 


83 ± 43 


62 ±46 




139 ± 38 


167 ± 38 


205 ± 39 




(154 ± 39) 


83 ± 17 


59 ± 17 


64± 16 


81 ± 15 


65 ± 17 


87± 15 


108 ± 16 


(106 ± 15) 


43 ±8 


53 ± 9 




74 ± 7 




(125 ± 6) 


(154 ± 8) 


(148 ±6) 


20 ± 19 


-2 ± 36 


36 ± 18 


15 ± 17 


72 ± 17 


46 ± 20 


70 ± 18 


(94 ± 17) 


64 ±33 






128 ± 27 


142 ± 28 






185 ± 28 


82 ± 16 


86 ± 9 


79 ± 13 


110 ± 7 


88 ±8 


105 ± 18 


126 ± 7 


(139 ± 7) 


69 ± 10 


74± 13 


89 ± 11 


91 ± 7 


97 ±9 


(124 ± 11) 


(132 ± 7) 


(151 ± 8) 



(14 ± 16) 

30 ± 18 
22 ± 15 
22 ± 13 
3 ± 12 
14 ±3 

(21 ± 8) 
11 ± 5 

31 ± 13 
(39 ± 17) 

(9 ± 13) 
(28 ± 10) 

52 ± 14 
29 ± 17 
46 ± 12 

32 ± 13 
-3 ± 13 

29 ± 8 
29 ± 10 
(47 ± 16) 

(49 ± 14) 
56 ± 11 

27± 10 
(58 ± 10) 
(39 ± 11) 



(45 ± 39) 
34± 17 
19 ±8 
3 ± 23 

40 ±8 
26 ± 11 



(-4 ± 16) 
44 ± 15 



27 ± 5 
61 ± 28 
(81 ± 8) 

34 ± 13 



39 ± 27 



25 ± 9 
42 ± 12 
72 ± 17 



68 ± 12 
52 ± 14 

86 ± 21 

31 ± 16 

11 ± 19 

45 ± 11 
25 ± 13 



29 ± 12 
56 ± 13 

33 ± 15 
37 ± 11 

18 ± 12 
24 ± 3 

107 ± 16 

34 ± 7 

19 ±4 
52 ± 12 

75 ± 16 
13 ± 13 

(50 ± 10) 
89 ± 8 
43 ± 17 
67 ±9 
51 ± 11 
61 ± 12 
54 ± 8 
45 ± 9 

76 ± 15 
88 ± 14 

73 ± 12 
68± 10 
54 ±21 

35 ±9 
88 ±9 
59 ± 8 

74 ± 11 

83 ± 15 

101 ± 39 
43 ± 16 

36 ± 7 
15 ± 18 
81 ± 30 

72 ± 7 
56 ±8 



25 



40 ±8 



52 



±4 



± 14 



101 ± 19 



66 ±38 
40 ± 15 



45 
(95 



±27 
±24) 



76 ±34 



Table 1 — Continued 





Ar I*" 














o r'" 














1048.220 


919.917 


922.200 


925.446 


916.815 


930.257 


919.658 


921.857 


924.950 


950.885 


976.448 


948.686 


971.738 


Star Name 


2.440 


-0.788 


-0.645 


-0.484 


-0.362 


-0.301 


-0.137 


-0.001 


0.155 


0.176 


0.509 


0.778 


1.052 






11.6 


15.3 


20.6 


25.0 


27.8 


34.5 


40.4 


46.9 


49.0 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


PG1605+072 


131 ± 15 






72 ±21 




99 ± 25 






136 ± 20 


163 ± 20 


229 ± 20 


(166 ± 22) 


(186 ± 20) 


PG1610+519 


122 ± 9 






81 ± 12 










150 ± 9 




256 ± 5 




(228 ± 7) 


PG2158+082 


117 ± 7 






100 ± 19 




97 ± 26 






121 ± 19 


175 ± 17 




(154 ± 22) 




PG2317+046 


68 ± 20 


49 ± 18 


29 ± 31 


61 ± 19 




50 ± 24 


28 ± 21 


48 ± 33 


71 ± 18 


135 ± 18 




132 ± 19 


(81 ± 19) 


Ton 102 


51 ± 13 


21 ± 15 


22 ± 16 


26 ± 16 




39 ± 16 


34 ± 16 


79 ± 16 


42 ± 17 


63 ± 16 


72 ± 15 






Ton S227 .... 


80 ± 25 


27± 10 


26 ± 11 


43 ± 10 




21 ± 20 


27 ± 11 


30 ± 12 


58 ± 9 




85 ± 11 


111 ± 12 


121 ± 9 


UV0904-02' 


69 ±4 


(30 ± 9) 




43 ±8 




59 ± 12 


52 ± 10 




67 ± 7 


65 ± 11 


75 ± 7 


84 ±8 


(89 ± 7) 



^Equivalent widths are given in mA. Values given in parentheses indicate lines that were not used in the linear fits. 

'^Numbers given below indicate (1) wavelengths in A, (2) line strengths in terms of log(/A) taken from Morton (2003), and (3) the equivalent width in mA for iV(0 I) = 

10^*^ cm~^ and 6 ^ 6 km s"^ 

*^Lines are arranged in order of increasing strength. 

^Name recognized by SIMBAD: 2MASS J20515997-4042465. 

^Central star of the planetary nebula NGC 6905. 

^Name recognized by SIMBAD: 2MASS J09070812-0306139. 



The sources of error mentioned above are 
straightforward to evaluate and would apply to 
just about any measure of an equivalent width. 
However, with the subdwarf stars, we must also 
contend with the confusion produced by stellar 
lines. We made no attempt to model such fea- 
tures, because in order to do so we would need 
to know the details of the stellar parameters for 
each star. Instead, we regarded the influence of 
stellar features as random sources of error in our 
line measurements. In order to estimate the am- 
plitude of such errors, which can vary markedly 
from one star to the next and can change with 
wavelength, we measured for each target the vari- 
ance of a large number of equivalent width mea- 
surements of imaginary, fake lines at wavelengths 
similar to the ones under study, but that were 
displaced away from known real interstellar lines, 
both atomic and molecular. This variance was 
then used as a guide for estimating the errors that 
should arise from stellar features. 

Figure [T] shows samples of spectra covering the 
relevant wavelength regions for two stars. These 
two cases illustrate strong differences in the degree 
of interference from stellar lines. The first exam- 
ple, AGK-l-81 266, has many stellar lines that can 
either add an apparent absorption to an interstel- 
lar line or distort the continuum level that is mea- 
sured on either side of the line. For this target, 
these effects dominate over other sources of error 
and create la uncertainties in W\ equal to 18.5 
and 14.2 mA for the Ar I and O I lines, respec- 
tively. This star also exhibits molecular hydrogen 
features of moderate strength, but here the lines 
are not strong enough to compromise the mea- 
surements of the atomic lines. A far more favor- 
able case for measuring interstellar features is pre- 
sented by the star UV0904— 02. Here, uncertain- 
ties produced by random stellar features should 
create errors of only 3.6 mA for the Ar I line and 
6.6 mA for the O I lines. 

All of the errors discussed in this section were 
combined in quadrature to synthesize the overall 
errors in the equivalent widths. Values for the 
equivalent widths of all lines and their associated 
uncertainties are listed for each of our targets in 
Table [T] The columns in this table are arranged 
in a sequence from the weakest to the strongest 
lines. 



4.2. Derivations of [Ar I/O I] Values and 
Their Uncertainties 

4-2.1. Reference Abundances 

Detailed discussions on various methods of 
measuring the protosolar and B-star abundances 
of Ar have been presented by Lodders (2008) 
and Lanz et al. (2008). We adopt a mean 
value for the recommended outcomes of the two, 
log(Ar/H) + 12 = 6.60 ± 0.10. (Since both de- 
terminations might be subject to common errors, 
the error of the mean is not reduced below the 
0.10 dex errors specified by each of them.) This 
value is higher than the solar photospheric value 
proposed by Asplund et al. (2009), and it remains 
so even after one applies a correction for gravi- 
tational settling of -1-0.07 dex (Lodders 2003) to 
obtain a protosolar value of log(Ar/H) -I- 12 = 
6.47 ± 0.13. For O, we take the solar photo- 
spheric value given by Asplund et al. (2009) 
and again apply a -1-0.07 dex settling correction 
to get log(0/H) -f 12 ^ 8.76 ± 0.05. This value 
agrees remarkably well with the measurement 
log(0/H) -I- 12 = 8.76 ± 0.03 obtained for B-stars 
by Przybilla et al. (2008). 

We must now consider the prospect that some 
of the Ar and O atoms are incorporated into solid 
form within or on the surfaces of dust grains, and 
this effect might be large enough to distort our 
findings on the differences in ionization. Unfortu- 
nately, we have no direct information about the 
depletion of gas-phase Ar, since ionization correc- 
tions (the object of the present study) can influ- 
ence the outcome. While in principle it would be 
beneficial if we could study iV(Ar I) along sight 
lines that penetrate dense media, where depletions 
are likely to dominate over ionization effects, this 
is not possible because the absorption lines are far 
too saturated (much more so than in the current 
study). 
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Fig. 1. Examples of FUSE spectra, where signals from all detector channels have bcicii combined, for the 
two stars AGK+81 266 (top group of panels) and UV0904— 02 (bottom group). Various interstellar features 
are indicated. The two panels showing the Ar I feature best illustrate the large contrasts in the strengths of 
unidentified stellar features that can interfere with the interstellar ones. 
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SJ98 presented a number of theoretical argu- 
ments that suggested Ar is not appreciably de- 
pleted in the low density ISM that we can observe. 
However, it would be good to confirm this out- 
look by some independent, experimental means. 
Fortunately, krypton is an element that can be 
observed in the ISM, and, like argon, is chemi- 
cally inert. It would be reasonable to expect that 
the capture of Kr onto interstellar dust grains, if 
it happens, would be similar to that of Ar. An 
advantage of studying Kr is that its interstellar 
features are weak (Cartledge et al. 2008), which 
means that they can be used to measure column 
densities over sight lines that have high values of 
nijltot) where element depletions should be gener- 
ally very strong. As with Ar, Kr has a photoioniza- 
tion cross section that is substantially larger than 
that of H (Sterling 2011). Thus, any simple mea- 
sure of the deficiencies of this element in the low 
density ISM could simply be a product of it being 
more easily photoionized than H. 

A way to overcome the confusion from the off- 
set produced by ionization is to compare differen- 
tial capture rates of elements onto grains as the 
conditions that favor grain formation change. For 
instance, Jenkins (2009) has determined that for 
sight lines with 7V(H°) > 10^^'^, where ionization 
corrections should be small, when 5 x 10^ O atoms 
are removed from the gas phase, 1(±1) atom of Kr 
vanishes. Since O is more abundant than Kr by 
a factor of 2.5 x 10^, any relative decrease in the 
abundance of Kr in the gas phase would be about 
half that of O (but the errors allow for this fac- 
tor to range from zero to being equal to that of 
oxygen). If we accept the idea that Ar depletes in 
the same manner as Kr, probably to within a fac- 
tor of \/mKr/m-Ar, and that in low density media 
O shows very low depletions (< 0.1 dex) (Jenkins 
2009), it is reasonable to adopt a gas-phase abun- 
dance ratio that is virtually the same as the proto- 
solar ones, log(Ar/O)0 = —2.16 ± 0.11. If indeed 
there is some mild depletion of O due to the for- 
mation of silicate dust grains, we may understate 
the strength of the ionization of Ar. 

4-2.2. Interpretation of Line Strengths 

We adopt the premise that the distribution of 
radial velocities of the neutral argon atoms is iden- 
tical to that of neutral oxygen (but this is not 
exactly correct; we will revisit this issue later). 



Figure [2] shows examples of some standard curve- 
of-growth plots for the O I lines appearing in the 
spectra of the same two targets that were featured 
in Fig.IH AGK+81 266 and UV0904-02. The for- 
mer of the two illustrates an average amount of 
line saturation for relevant features, while the lat- 
ter represents an extreme case of saturation caused 
by a low overall dispersion of radial velocities. In 
principle, we could have derived values for A^(0 I) 
from the best-fit curves of growth shown by the 
dashed curves in the figure panels and then assume 
that 7V(Ar 1) follows from the equivalent width of 
the one available line at 1048.220 A assuming the 
same velocity dispersion parameter b as that found 
for O 1. Instead, we used a much simpler approach 
that sidesteps the goal of deriving explicit values 
of N and b (whose errors are strongly correlated) 
and proceeds directly to an answer for just the ra- 
tio of the two column densities. An advantage here 
is that we can make a straightforward analytical 
determination of the uncertainty of the outcome 
based on the errors of some linear fitting coeffi- 
cients. 

The comparison of Ar I to O 1 is based on the 
following principle. If one could imagine the exis- 
tence of a hypothetical O I line with a transition 
strength log(/A)o i that is just right to produce 
a value of W\ / A that exactly matches that of the 
Ar I line, one could then derive the deficiency of 
Ar I with respect to its expectation based on O I, 
which we denote in logarithmic form as [Ar I/O I]. 
This quantity yields the logarithm of the ratio of 
the two neutral fractions relative to the solar abun- 
dance ratio and is given by the relation 

[Arl/OI] = log(/A)o I - log(/A)Ar I - log(Ar/O)0 
= log(/A)o i-0.28±0.11, (1) 

where log(/A)Ar i = 2.440 ±0.004 (Morton 2003). 
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Fig. 2. — Examples of how the values of [Ar I/O I] are derived for the two stars chosen for Fig. [T] The 
equality shown in each x-axis label applies to the horizontal projection of the Ar I line strength onto the O I 
curve of growth, as expressed in Eq. [TJ We performed weighted least-squares linear fits for the logarithms 
of Wx/X for lines of O I that are most influential in establishing the trends (solid lines) for a comparison 
with the measurements of \ogW\/X of the Ar I line. Values for cq and ci given in Eq. [2] that define these 
trends are shown in the boxes. The best-fit curves of growth together with their b values are also shown, 
but they are not used in the derivations. O I lines that had strengths that were well above or below those 
in the important portion of curve of growth were not included in the derivation of the best fits, and they 
are indicated here as "(not used)." These unused lines are identified with parentheses around the equivalent 
width values listed in Table [1] 
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To determine the strength of the hypothetical 
O 1 line, we perform a weighted least-squares Un- 
ear fit for log(W^A/A) vs. the quantity on the right- 
hand side of Eq. [T] (the abscissa for each plot in 
Fig. [5]) for an appropriate selection of O I lines. 
[As with the Ar I line the /-values of the O I tran- 
sitions are from Morton (2003).] The lines that 
are chosen for this fit are ones that are situated 
not too far from the horizontal projection (shown 
by dotted lines in the figure) of log(W^A/A) for the 
single line of Ar I onto the trend for the O I lines. 
With this restricted fit, we define a simple relation- 
ship in log(WA/A) that is a good approximation to 
a relevant portion of the curve of growth. 

In the two panels of Fig. [21 these best-fit linear 
trends are shown by the straight solid lines. They 
depict the relation between y — log Wx/\ and x = 
log(/A) — 0.28 according to the equation 

y ^ cq + ci{x - xq), (2) 

where 

x,^Y.^,a{y)f/Y,a{y)-^ (3) 

i i 

represents a zero reference point that produces a 
vanishing covariance for the errors in the fitting 
coefficients cq and ci. A nominal value on the 
x axis for the projection of j/Ar i onto the linear 
relation is given by 



In the fraction part of this equation, the numer- 
ator and denominator have errors (cr(yAr i)^ + 
(t(co)^)°'^ and cr(ci), respectively. A conventional 
approach for deriving the error of the quotient is 
to add in quadrature the relative errors of the two 
terms, yielding the relative error of the quotient. 
However, this scheme breaks down when the er- 
ror in the denominator is not very much less than 
the denominator itself. A more robust way to de- 
rive the error of a quotient has been developed 
by Geary (1930); for a concise description of this 
method see Appendix A of Jenkins (2009). We use 
this method here; it is effective as long as there is 
little chance that the denominator minus its error 
could become very close to zero or be negative, 
i.e., denom./cr(denom.) > 3. 

The dotted lines in Figure [5] show schematically 
how the best fit values and the error ranges for 



[Ar I/O I] are derived. Note that the horizontal 
and vertical segments for the error limits do not 
exactly intersect the best linear trend for the O I 
lines because the error analysis allows for the un- 
certainty for the location of this line. (But we 
point out that the intersections for the worst pos- 
sible error in one direction for yxr i do not occur 
at the locations for the worst possible deviations 
in the opposite direction for the trend line.) Also, 
the final errors for [Ar I/O I] — XAr i are not sym- 
metrical about the best values. On average, the 
upward error bounds are about 75% as large as 
the negative ones. 

The difference in atomic weights of Ar and 
O will cause the thermal contributions to the 
Doppler broadenings of these two elements to dif- 
fer from each other. The impact of this effect on 
our results for the column density ratios is small 
however. For instance, if we consider that we are 
viewing absorption lines arising from the warm, 
neutral medium (WNM) and there were no bulk 
motions of the gas, the line broadening parame- 
ters fothoim. caused by thermal Doppler broadening 
for T = 7000 K would be 2.7 and 1.7km s"^ for 
O 1 and Ar I, respectively. For a typical observa- 
tion, such as the one for AGK-f 81 266 illustrated 
in the left-hand panel of Fig. [2l we find that the 
observed curve of growth for the O 1 lines, yield- 
ing an apparent 6obs. = 9.9km s^^, indicates that 
kinematic effects arising from turbulent motions 
(or multiple velocity components) should be the 
most important contribution to the broadening 
since 6turb. = V^^Z^WW^. = 9-53 km s-^ 
is only slightly smaller than 6obs.- The curve of 
growth that characterizes the saturation of the 
Ar I line would conform to a slightly lower ve- 
locity dispersion parameter compared to that ob- 
served for^O I, 6ob.. = x/^turb. + ^(Ar I)Lrm. = 
9.68 km s ^, because the higher atomic weight of 
Ar causes the thermal contribution to be smaller. 
For the observed equivalent width of the Ar I line, 
the error in the ratio of column densities caused 
by our assumption that the values of 6obs. of the 
two elements are identical will create an underesti- 
mate of [Ar I/O I] = — 0.012 dex. A similar calcu- 
lation for the more extreme line saturation exhib- 
ited by UV0904— 02 shown in the right-hand panel 
of Fig. [2] indicates that the perceived outcome for 
[Ar I/O I] could be too low by —0.12 dex. Cases 
showing this much saturation are rare for our col- 
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lection of sight lines. 

If the kinematic line broadening is not an ap- 
proximately Gaussian form that one would ex- 
pect from pure turbulent broadening, but instead 
results from distinct and well-separated narrow 
components, the errors in [Ar I/O I] could be 
larger than those evaluated above. However, the 
recordings of Ar I lines for 9 different stars made 
by IMAPS at a resolution of 4 km s~^ that were 
shown by SJ98 reveal profiles that, while not ex- 
actly Gaussian, are nevertheless generally smooth 
and devoid of any narrow spikes that are isolated 
from each other. Thus, the numerical estimates 
presented in the above paragraph should be rea- 
sonably accurate. 

One might question whether or not errors in 
the adopted /-values could cause global system- 
atic errors in the evaluations of [Ar I/O I]. While 
Morton (2003) listed a very small uncertainty for 
the /-value of the Ar I line at 1048.220 A, he did 
not specify errors for the O I lines. Nevertheless, 
empirical evidence from high quality FUSE ob- 
servations of WD stars by Hebrard et al. (2002), 
Sonneborn et al. (2002) and Oliveira et al. (2003) 
all showed curves of growth that are remarkably 
well behaved for the same lines that are used in 
the present study. While one could still pose the 
objection that all of the O I lines could collec- 
tively have a systematic error of a certain magni- 
tude, and yet could still yield acceptable curves of 
growth, this seems unlikely: the values of [O I/H I] 
derived by Sonneborn et al. (2002) and Oliveira 
et al. (2003) are generally consistent with those 
found elsewhere in the ISM based on measure- 
ments of the intersystem O I line at 1355.6 A 
(Jenkins 2009). [Hebrard et al. (2002) did not 
attempt to measure iV(H I).] 

4-2.3. Outcomes 

Table [5] shows the outcomes of our analysis for 
all of the targets in the survey, along with the ap- 
plicable Galactic coordinates and apparent mag- 
nitudes of the stars. The last two columns present 
some cautions in numerical form. First, Column 
(9) lists for the fraction part of Eq. 2] the denom- 
inator divided by its error. If this number is less 
than about 3, the upper limit for [Ar I/O I] should 
not be trusted. Second, Column (10) lists the 
probability of obtaining a worse fit of the O I lines 
to the linear trend (i.e., a higher value for x^), 



given the errors that we derived. A plot of the 
frequency of all of these numbers shows a distribu- 
tion that is consistent with a uniform distribution 
between and 1, which indicates that our error 
estimates for W\ of the O I lines are neither too 
conservative nor too generous. Thus, any individ- 
ual case where this probability is low should not 
be considered as a real anomaly. 
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Fig. 4. — Locations of the stars and their values of 
[Ar I/O I] shown on an AitofF projection of Galac- 
tic coordinates. The sizes of the black circles indi- 
cate the values of [Ar I/O I] according to the key 
shown at the bottom. 

Figure [3] shows a graphic representation of all 
of the results and their uncertainties. They were 
ranked and then arranged in order of the best val- 
ues of [Ar I/O I] to make it easier to see the dis- 
persion of results and also to show that the more 
extreme deviations often represent cases where the 
errors are somewhat larger than normal. 
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Table 2 
Results for [Ar I/O I] 
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LBS 1362 


273.67 


6.19 
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12.50 


-0, 


.60 


-0.53 


-0.48 


7.05 


0, 


.65 


MCT 2005-5112 


347.71 


-32.68 


15.20 




-0, 


.63 


-0.47 


-0.36 


3.82 


0, 


.46 


MCT 2048-4504 


1.10 


-39.51 


14.85 


15.20 


-0, 


.52 


-0.28 


-0.09 


4.97 


0, 


.97 


NGC6905 star 


61.49 


-9.57 


16.30 


14.50 


-0, 


.46 


-0.26 


-0.12 


2.56 


0, 


.75 


PG0919+272 


200.46 


43.89 


12.27 


12.69 


-2, 


.24 


-0.89 


-0.16 


3.83 


0, 


.64 


PG0952+519 


164.07 


49.00 


12.47 


12.80 


-0, 


.60 


-0.51 


-0.44 


5.11 


0, 


.90 


PG1032+406 


178.88 


59.01 


10.80 




-3, 


.79 


-1.84 


-0.41 


2.71 


0, 


.30 


PG1051+501 


159.61 


58.12 


14.59 




-0, 


.01 


0.28 


0.59 


2.77 


0, 


.61 


PG1230+068 


290.42 


68.85 


12.25 




-0, 


.57 


-0.23 


0.06 


7.58 


0, 


.03 


PG1544+488 


77.54 


50.13 


12.80 




-0, 


.48 


-0.40 


-0.34 


5.41 


0, 


.56 


PG1605+072 


18.99 


39.33 


13.01 


12.84 


-0, 


.52 


-0.36 


-0.24 


5.31 


0, 


.88 


PG1610+519 


80.50 


45.31 


13.73 




-0, 


.53 


-0.44 


-0.36 


10.93 


0, 


.22 


PG2158+082 


67.58 


-35.48 


12.66 




-0, 


.91 


-0.64 


-0.49 


2.78 


0, 


.16 


PG2317+046 


84.84 


-51.10 




12.87 


-1, 


.42 


-0.79 


-0.37 


4.20 


0, 


.04 


Ton 102 


127.05 


65.78 


13.29 


13.54 


-0, 


.89 


-0.45 


-0.15 


2.96 


0, 


.31 


Ton S227 


201.39 


-77.81 


11.60 


11.90 


-0, 


.34 


0.12 


0.44 


9.06 


0, 


.36 


UV0904-02 


232.98 


28.12 


11.64 


11.96 


-0, 


.39 


-0.19 


-0.03 


3.78 


0, 


.93 



''[Ar I/O I] = — (Ar/O)0 — T^i^e range of possible values shown here do not include an overall systematic error that arises 
from the uncertainties in the constants that are included in Eq. [T] which amounts to a combined error of 0.11 dex. 

'^Coordinates and apparent magnitudes were supplied by the SIMBAD database. 

'^The relevance of this quantity is discussed in the text that follows Eq. |4] If the listed value is less than about 3, the positive 
value of the error quotient may be misleading. 
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Figure |4] shows the locations of the targets in 
the sky and their respective vahies of [Ar I/O I]. 
While no obvious regional trends seem to be evi- 
dent, we can perform a test to determine whether 
or not the variability of the outcomes exceeds what 
we would have expected from our errors (assuming 
that they are correct). We do this by computing a 
value for , where the choice for the error of each 
case depends on whether a test value is above or 
below the measurement outcome. This procedure 
properly takes into account the asymmetries of the 
errors. Adopting this method, we find that a min- 
imum of 66.9 for the 44 measurements occurs 
at a value [Ar I/O I] = —0.438. This minimum 
for the with 43 degrees of freedom is greater 
than what we would have expected from our er- 
rors alone (the probability of obtaining by chance 
a higher value for x^ is only 1%). If we were to 
propose that the real variability across the sky is 
a = 0.11 dex, and we add this value in quadrature 
to all of the experimental errors, the minimum x^ 
drops to 42.4, which makes the probability of a 
worse fit equal to 50%. The location of this new 
minimum is at [Ar I/O I] = —0.427. 

The discussion above has considered only the 
random fluctuations arising from either the mea- 
surements or the true variability in [Ar I/O I] in 
different sight lines. We must not lose sight of the 
fact that there can be an overall systematic error 
of 0.11 dex for the entire collection. This global 
error arises from the uncertainties of the value for 
log(Ar/O)0 that went into Eq.[T] It is much larger 
than the error in the weighted mean value for all 
of the measurements. 

5. ARE THE STARS BEYOND THE 
BOUNDARY OF THE LOCAL BUB- 
BLE? 

As discussed in Section [21 our objective is to 
sample interstellar material that is beyond the 
edge of the Local Bubble, enough so that our mea- 
surements are not heavily influenced by the very 
low density gas within this cavity. In principle, we 
could compare the three-dimensional locations of 
the target stars with maps that outline the bound- 
ary of the Local Bubble (Vergely et al. 2010; 
Welsh et al. 2010; Reis et al. 2011) to indicate 
whether or not we are primarily sampling gas in 
the surrounding denser medium. However, the dis- 



tances to our targets are uncertain, which makes 
this approach unworkable. Instead, we adopt a 
definition of the boundary proposed by Sfeir et al. 
(1999) (and one that was also used by Lehner et 
al. (2003)), who linked its location to the sud- 
den onset of Na I D-line absorption that crossed 
a threshold VFa(D2) = 20 mA. This threshold is 
equivalent to A^(H I) « 2 x lO^^cm-^ (Ferlet et 
al. 1985), which in turn corresponds to A'^(0 I) 
that is slightly greater than lO^^cm"^. Lehner et 
al. (2003) found that a typical velocity disper- 
sion parameter 6(0 I) = 6 km s~^ occurred inside 
the Local Bubble. Using these two parameters for 
O I, {N, b) = (10^^cm~^, 6km s~^), we can com- 
pute the equivalent widths for all of the O I lines 
when the boundary is crossed. The third row of 
numbers in the column headings in Table [1] shows 
the values of W\ (in mA) for all but the strongest 
three lines. By comparing these values with the 
entries that show our measurements, particularly 
the weaker but securely measured ones, we can as- 
certain that our targets are beyond the edge of the 
Local Bubble. 

6. INTERPRETATION: FUNDAMEN- 
TAL PROCESSES AND EQUATIONS 

In this section, we address the basic physical 
processes that relate our findings on [Ar I/O I] to 
the ionization balance of the gas and the resulting 
degree of partial ionization. Our discussion about 
the means of ionizing the atoms will focus mainly 
on the primary ionization by photons, along with 
the effects of collisional ionizations caused by sec- 
ondary electrons that originate from these primary 
ionizations. These two processes are the most im- 
portant sources of ionization, and they represent 
one side of the balance between recombinations 
with free electrons and charge exchange reactions 
between various constituents of the medium. For 
completeness, we will also cover other means of 
ionizing atoms and creating free electrons, such 
as cosmic ray ionizations, the ionizations of in- 
ner shell electrons of heavy elements by X-rays, 
the creation of ionizing photons when helium ions 
recombine, and the nearly complete ionization of 
many elements that have first ionization potentials 
below that of hydrogen. 
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6.1. Direct and Secondary Ionizations of 
H and Ar 

6.1.1. Photoionization 

The primary photoionization cross sections for 
neutral Ar are larger than those for H by about one 
order of magnitude at low energies, and the ratio 
increases substantially at higher energies. Various 
secondary ionizing processes initiated by the pri- 
mary photoionizations of H and He likewise have 
a stronger effect on Ar than on H. This contrast 
in ionization rates is the fundamental tool that we 
use in the interpretation of the Ar data to quantify 
the photoionization of H and the subsequent cre- 
ation of free electrons. For the collective effect of 
all of these ionization channels, we can construct a 
simple formalism based on arguments created by 
SJ98. They defined a quantity based on ioniza- 
tion rates F and recombination coefficients a for 
the two elements, 



F(Ar°)Q(H°,T) 
F(H")a(Ar",T) 



(5) 



SJ98 considered only primary photoionizations of 
these two elements. Going beyond their develop- 
ment, we construct a more comprehensive picture 
by considering some refinements in the calcula- 
tions of F for both H° and Ar°. 

First, we start with the primary photoioniza- 
tion rates Fp = / F{E)a{E)dE, where F{E) is 
the ambient photon flux as a function of energy 
E and a{E) is the photoionization cross section 
for either (Spitzer 1978, pp. 105-106) or Ar° 
(Marr & West 1976). Next, we include secondary 
ionizations with rates F^ that are created by the 
collisions from energetic electrons that are liber- 
ated by the primary photoionizations of H and 
He. Added to this are the effects from photons 
with energies above about 300 eV, which can inter- 
act with the abundant heavy elements in the ISM 
to produce additional energetic electrons that will 
ionize H and Ar with a rate that we identify as F^/ . 
These electrons arise from the primary ionizations 
of the inner electronic shells, and they are supple- 
mented by one or more additional electrons from 
the Auger process. Finally, we must acknowledge 
that recombinations of singly- and doubly-ionized 
He ions with electrons create additional photons 
when de-excitation occurs in the lower stage of 
ionization, either He*^ or He^. The importance of 



not overlooking secondary electrons and recombi- 
nation photons from the He ionizations is under- 
scored by the fact that while the abundance of He 
is only 1/10 that of hydrogen, its primary ioniza- 
tion cross section is 6 to 100 times that of H over 
the energy range 25 to 4000 eV. Most of the re- 
combination photons are capable of ionizing both 
H and Ar, and they supplement the other sources 
of ionization with rates Fji/e" and THe+ ■ In short, 
we consider that the total ionization rates for the 
two elements in Eq. [5] each consist of five contri- 
butions. 



F-F, 



He+ 



(6) 



The details of how we compute F^, Fs', F//e+ and 
FjLf gO are discussed in Appendices |X] and |B1 

The quantity Pai- defined in Eq. [5] provides a 
means for evaluating how large the neutral frac- 
tion of Ar should be relative to that of H according 
to the formula 



[Ar I/H I] = log 



l + n(H+)/n(HO) 



1 + PArn(H+)/n(HO) 



(7) 



For a more accurate formulation of [Ar I/HI] that 
will be developed later in Section 16. 4[ we will in- 
troduce a more refined parameter P^^. that will be 
based on a calculation that is more elaborate than 
the one shown in Eq. [5l This new parameter will 
be substituted for Pat in Eq. [T] Under most con- 
ditions, the differences between P^^, and Pai- are 
small. 

Figure [5] shows for both H and Ar the mo- 
noenergetic cross sections for primary ionization 
(dashed lines), the effective enhancements arising 
from the secondary ionizations Tg (dash-dot lines) 
for Xe — 0.05, and the additional effects from F^', 
^He+ and TheO i all of which give total ionization 
rates shown by the solid lines. For H, secondary 
ionizations outweigh the primary ones for photon 
energies above about 100 eV. By contrast, we find 
that for Ar the enhancement at these energies is 
small compared to its much higher primary cross 
section. A more narrowly defined version of Paf, 
which we denote as Px-c{E), applies to an irradi- 
ation of the gas by photons of a given energy E, 
rather than from a combination of fluxes over a 
broad energy range. The definition of Pat{E) is 
illustrated in the upper panel of the figure, and its 
dependence on E is shown in the bottom panel. 
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Fig. 5. — Upper panel: Various combinations of effective photoionization cross sections for neutral hydrogen 
and argon as a function of energy E. The dashed Unes depict the absorption cross sections for the primary 
ionization rates Tp that would apply to monoenergetic radiation at the energy E shown on the a;-axis. The 
effective cross sections that arise by including the additional effects from ionizations by secondary electrons 
produced by H and He (related to Tg) are shown by the dot-dash lines. These secondary ionization efficiencies 
were calculated according to the principles outlined in Section lA.ll of the Appendix for an electron fraction 
Xe — 0.05.; see Eqs. lAll and IA21 These two ionization processes can be further enhanced by photons 
that come from recombinations of He ions with free electrons (Appendix |B| and electrons from inner shell 
ionizations of heavy elements (Appendix IA.2I) , leading to total effective cross sections depicted by the solid 
lines. Across the top of the plot are markers showing the energies at which half-intensity penetration depths 
for the photons occur for various hydrogen column densities, assuming the gases have a solar abundance 
ratio for various elements. The quantity PAr{E) defined in Eq. [5] for a specific energy E is approximately 
equal to the ratio of the total effective photoionizatio't? cross sections that give rise to r(H°, Ar°) defined in 
Eq. ini since the recombination rates for the two elements are about the same for all temperatures. Bottom 
panel: A plot of Pat{E) as a function of energy. 



6.1.2. Other Ionization Mechanisms 

As discussed earlier in Section [U neutral atoms 
in the ISM are subjected to ionization by cosmic 
ray particles. These ionizations, with a rate CcRj 
add to the effects of primary and secondary ion- 
izations from photons discussed above. To ob- 
tain proper values of Ccr that apply to the dif- 
fuse medium, we must apply corrections to the 
values of Ccr that were measured for the more 
dense media that have appreciable concentrations 
of molecules. The details of this computation are 
discussed in Appendix [Cl 

In the WNM, the main source of free electrons 
is from the ionization of H and He. However, a 
small number of additional electrons arise from 
other atoms that can be ionized by starlight pho- 
tons less energetic than the ionization potential of 
H. For the combined effect from these elements, we 
adopt an estimate n(M+) equal to 1.5 x lO"** times 
the density of hydrogen based on the assumption 
that the gas we are viewing is in a regime where 
the element depletions are relatively modest. 

6.2. Recombination 

The radiative recombination coefficients for free 
electrons and ions to create neutral hydrogen and 
singly-ionized helium, a(H*',T) and a(He^,T), 
are taken from Spitzer (1978, pp. 105-107). For 
a(H°), we excluded recombinations to the low- 
est electronic level, since they generate Lyman 
limit photons that can re-ionize hydrogen atoms 
over a short distance scale. Recombination co- 
efficients for He° were taken from Aldrovandi & 
Pequignot (1973) and those for Ar° from Shull & 
Van Steenberg (1982). The results for Ar° given 
by Shull & Van Steenberg (1982) agree well with 
the radiative recombination coefficients listed by 
Aldrovandi & Pequignot (1974). The minimum 
temperature where dielectronic recombination for 
Ar° becomes important is 2.5 x 10** K (Aldrovandi 
& Pequignot 1974), which is above temperatures 
that we will consider, hence we can ignore this 
process. 

In addition to recombining with a free elec- 
tron, an ion can also be neutralized by colliding 
with a dust grain and removing an electron from 
it (Snow 1975; Draine & Sutin 1987; Lepp et al. 
1988; Weingartner & Draine 2001a). The opera- 
tion of this effect on protons is important for regu- 



lating the fraction of free electrons in the cold neu- 
tral medium (CNM), but it is of lesser significance 
for the WNM [cf. Figs. 16.1 and 16.2 of Draine 
(2011)], which probably dominates the sight lines 
in our study. The rate constant ag for this pro- 
cess is normalized to the local hydrogen density 
'T-(IItot) = ri(H°) -I- n(H+), and it depends on sev- 
eral physical parameters that influence the charge 
on the grains, such as the electron density n(e), 
the rate of photoelectric emission that is driven 
by the intensity G of the local radiation field be- 
tween 6 and 13.6 eV, and the temperature T. For 
our equilibrium equations in Section 16.41 we have 
adopted parametric fits for ag(H°, n(e), G, T) and 
ag(He°, n(e), G, T) from Weingartner & Draine 
(2001a). They do not supply fit coefficients for 
ag(Ar°, n(e), G, T), but since the ionization poten- 
tial of Ar° is close to that of H'^ , it is reasonable to 
adopt the hydrogen rate coefficient and divide it 
by the square-root of the atomic weight (40) of Ar. 
Throughout our analysis, we set G = 1.13, which 
is the value recommended for the general ISM by 
Weingartner & Draine (2001a). 

6.3. Charge Exchange 

In addition to the recombination processes 
mentioned in the previous section, charge ex- 
change reactions with neutral hydrogen can also 
lower the ionization state of an atom. With refer- 
ence to such reactions for element X, X+ -I- H° ^' 
XO -I- H+ and X++ + H" ^ X+ + H+, we adopt 
the notation G'(X+,r) and G'(X++,r) for the 
respective rate constants. For our calculations 
of G'(He+,T), G'(He++,r), G'(Ar+,r), and 
G'(Ar"'"+,T), we adopted Kingdon & Ferland's 
(1996) fits to the calculations from various sources 
(see their Tabic 1 for the coefficients and refer- 
ences). 

Since the ionization potential of neutral oxygen 
is close to that of hydrogen, the rate constant for 
charge exchange of these two species is large and 
reverse endothermic reaction is not negligible, ex- 
cept at very low temperatures. The rate constant 
G(0+,r) for the reaction 0° -h H+ ^> 0+ + H'' 
can be obtained from G'(0"'",T) by the principle 
of detailed balancing. However, in doing so, one 
must treat the three fine-structure levels of the 
ground state of 0° separately, since their energy 
separations are comparable to the differences in 
the ionization potentials of H and O. The large 
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rate constants in both directions (Stancil et al. 
1999) assures that the ionization fraction of O is 
locked very close to a value 8/9 times that of H 
for T > 10"^ K. This is a key principle that allows 
us to use O as a substitute for H in the present 
study or Ar vs. H fractional ionizations. 

For completeness, we should also consider 
charge exchange reactions between Ar and He, 
even though the abundance of He is much lower 
than that of H. The charge transfer recombina- 
tion reaction Ar"*"^ + He" At^ + He''" has a 
rate constant D'{Ai++,T) = 1.3 x lO-^Ocm^s-i 
according to Butler & Dalgarno (1980), which 
they claim to be constant for T > 10^ K. The 
charge transfer ionization reaction with He+, i.e., 
Ar'' + He^ Ar^ + He'', has a rate constant 
D{Ai+,T) < 10~^^cm3s~^ according to Albrit- 
ton (1978), so we will ignore this process in the 
equation for the ionization balance for Ar (Eq. |5] 
below). 

6.4. Equilibrium Equations 

In a medium where both hydrogen and helium 
are partly ionized, the densities of an element X in 
its 3 lowest levels of ionization X°, X"*" and X++ 
are governed by the equilibrium equationtH 



r(X°) + Ccr(X°) + C(X+, T)n(H+)J n(X'') 
a(X°, r)n(e) + C'iX+ ,T)n{R°) 
ag(XO, n(e), G, T)n(Htot)l n(X+) (8) 



and 



state with free electrons as a function of temper- 
ature T. The simultaneous solution to these two 
equations yields the fractional abundances in the 
three ionization levels 



/o(X,T) 



n(X") 



n(XO) +n(X+) +n(X++) 



= 1 





r(xO) 


+ CcR(X°) + C(X+,T)n(H+)' 




r(x+) + y 




a(XO,T)n(e) 


+ C'(X+, T)n(HO) + agiXO, n{e),G, r)n(Htot) 


Y 



(10a) 



with 



Y = a(X+,r)n(e)+C'(X++,r)n(HO)+i?'(X++,T)n(He'') 

(10b) 



/++(x,r)E 

and 

/+(X,T) = 



i(X++) 



^ i-./o(x,r) 

n(X") + n(X+) + n(X++) 1 + Y/T{X+) ' 

(11) 



n{X+) 



n{X") + n{X+) + n{X++) 



l-/o(X,T)-/++(X,r). 
(12) 



6.5. Electron Density and the Ionization 
Fractions of H and He 

Before the ionization fractions of Ar can be de- 
rived, we must determine not only the ionization 
balance of hydrogen, but also that for helium. We 
do this by solving Eqs. [TOl4T2l (substituting He for 
X and eliminating the D' term) along with the 
equation for the hydrogen ionization balance. 



r(X+)n(X+)= \a{X+,T)n{e) + C"(X++,T)n{B°) 



r(H")+CcR(H") + Z 



+ D'(X++,r)n(He°) n(X++) 

where r(X^) is the photoionization rate of element 
X in its ionization state y (neutral, -t-, or H — h) and 
a(X^,T) is the recombination rate of the y + 1 



(9) 



n(HO) a(HO, T)n(e) + agiX^, n(e), G, r)n(Htot) - Z ' 



(13a) 



with 



Z = 0.1n(HO)[G'(He+,T)/+(He,T)+G'(He++,r)/++(He,r)] , 

(13b) 

and the constraints 



^The development here follows that given by Eqs. 12—19 of 
SJ98, except that we have added the He charge exchange re- 
combination reaction as an additional channel for reducing 
the ionization of element X in its doubly charged form. We 
have also added cosmic ray ionizations and have implicitly 
included the various kinds of secondary photoionizations in 
the definition of F, as indicated in Eq.[6] We have corrected 
Eos. Island 1101 here to include a missing n(H+) term, which 
was a typographical oversight in the equations of SJ98. 



n(He")-fn(He+)-fn(He++) = 0.1n(H°) 



z(H+) 



and 



i(e) = n(H+) + n(He+) + 2n(He 



n(HO) 
(14) 



n(M+) . 
(15) 
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Since the hydrogen ionization balance depends 
on n(e) which in turn is influenced by the ioniza- 
tion fractions of He (which are also influenced by 
n(e)), we must solve Egs. ll3allT5] iterativelv to ob- 
tain a solution. We found that these equations 
converged very well if we kept n(Htot) and the 
ionization rate of H pegged to a certain valuej^l 
Starting with a zero helium ionization rate, the 
iterations progressed slowly to successively higher 
rates until the final, correct value was reached and 
the ionization fractions had stabilized. 

After obtaining the final results for the coupled 
hydrogen and helium ionization balances, we can 
solve for /o(Ar) using Eq.[Tni This result can then 
be used to derive a more accurate value for Pai-, 

which can be substituted for Pai- in Eq. [7] to ob- 
tain a solution for [Ar I/H I] that makes use of 
all of the physical processes that were discussed in 
Sections 

7. OUTCOME FROM KNOWN SOURCES 
OF PHOTOIONIZATION 

7.1. External Radiation 

Over many decades, the diffuse, soft X-ray 
background has been measured by a large num- 
ber of different experiments [for a review, see Mc- 
Cammon & Sanders (1990)]. Most of the emis- 
sion below 1 keV arises from hot (T > 10^ K) 
gas in the Galactic disk and halo, with radiation 
from extragalactic sources dominating at higher 
energies (Chen et al. 1997; Miyaji et al. 1998; 
Moretti et al. 2009). Much of the literature on 
the diffuse radiation shows a distinction between 
contributions from a local component with little 
foreground absorption and more distant emissions 
with varying levels of absorption. The local back- 
ground was once identified as having originated 
from hot gas in the Local Bubble (Sanders et al. 
1977; Hayakawa et al. 1978; Fried et al. 1980), 
but in recent years it has been recognized to be 
strongly contaminated, or completely dominated, 
by X-rays arising from charge exchange produced 
by the interaction of the solar wind with incoming 



Fixing n(e) instead of n(Htot) produces unstable, oscillat- 
ing solutions under certain circumstances. 



interstellar atoms (Cravens 2000; Lallement 2004; 
Pepino et al. 2004; Koutroumpa et al. 2006, 2007, 
2009; Peek et al. 2011; Crowder et al. 2012). For 
this reason, we ignore the weakly absorbed, nearly 
isotropic portion of the X-ray background and fo- 
cus our attention to the component that exhibits 
a pattern in the sky that clearly shows absorption 
by gas in the Galaxy. 

Kuntz & Snowden (2000) have performed a 
detailed investigation of the nonlocal compo- 
nent, which they call the transabsorption emis- 
sion (TAE). They describe the strength and spec- 
tral character of the TAE in terms of emissions 
from optically thin plasmas at two different tem- 
peratures. They define a soft component that 
has a mean intensity over the sky / = 2.6 x 
10~®erg cm~^s~"'^sr~-'^ over the interval 0.1 < E < 
2 keV and a spectrum consistent with the emission 
from a plasma at a temperature T — 10^"^ K, and 
this flux is accompanied by a hard component with 
/ = 8.5 X lO^^erg cm~^s~^sr~^ over the same en- 
ergy interval with T = 10'^ '"'K. To translate the 
sum of these two components into a distribution of 
the photon flux as a function of energy, F(E), we 
calculate synthetic flux representations using the 
CHIANTI database and software (Version 6.0) 
(Dere et al. 1997, 2009), after normalizing the 
emission measures to give the intensities stated 
above (we find that EM = IQi^-^'^ and W^^-^^cm-^ 
for the soft and hard components, respectively). 
We supplement the TAE result with an underly- 
ing power-law extragalactic emission of the form 
10.5 phot cm"2s-isr-ikeV~^i;(keV)-i ''6 (Chen 
et al. 1997). 

The ISM is opaque to X-rays at the lowest 
energies. The energies at which half of the X- 
rays are absorbed for various column densities 
are shown in the top portion of Figure [SI which 
were derived from the calculations of Wilms et al. 
(2000). At energies of around 100 eV where the 
ISM is neither completely opaque or transparent 
for iV(II°) « few X lO^^cm"^, uncertainties in the 
layout of emitting and absorbing regions make it 
difficult to calculate with much precision how far 
the X-rays can penetrate the typical gas volumes 
that were sampled in our survey of Ar I and O I. 
Thus, rather than implement an elaborate atten- 
uation function that would be difficult to explain 
(and perhaps not especially correct at our level of 
understanding), we apply a simplification that all 
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of the X-rays are transmitted above some thresh- 
old energy and none below it. The threshold that 
we adopted was 90 eV, on the assumption that in 
some directions the gas can view the unattenu- 
ated X-ray sky through 7V(H°) slightly less than 
IQi^cm-^. 

The upper panel of Figure [5] shows our synthesis 
of the sum of the extragalactic power-law emission 
and TAB synthesis described above. For the pur- 
pose of calculating F for various elements, we make 
use of only the flux depicted by the dark trace in 
the figure, i.e., that which starts at the cutoff en- 
ergy (90 eV) and ends at an energy beyond which 
no appreciable additional ionization occurs. 

7.2. Internal Radiation Sources 

Embedded within the ISM are sources of EUV 
and X-ray radiation that can make additional con- 
tributions to F. We can make estimates for their 
average space densities and the character of their 
emissions, but one uncertainty that remains is how 
well the ensuing photoionizations are dispersed 
throughout the ambient gas. At one extreme 
representing minimum dispersal, we envision the 
classical Stronigren Spheres that surround sources 
that are not moving rapidly and that emit most 
of their photons with barely enough energy to ion- 
ize hydrogen. These photons have a short mean 
free path in a neutral medium. Under these cir- 
cumstances, the zone of influence of the source is 
sharply bounded, and the resulting ionization is 
nearly total inside the region and zero outside it. 
At the other extreme, one can imagine that the 
photons, ones that have relatively high energy, can 
travel over a significant fraction of the inter-source 
distances before they are absorbed. In addition, 
the sources themselves could move rapidly enough 
that they never have a chance to establish a sta- 
ble condition of ionization equilibrium. (This issue 
will be investigated quantitatively in Section [731) 
These conditions could lead to the ionization be- 
ing more evenly distributed and not necessarily 
complete. If the sources have large enough veloc- 
ities, we can even imagine a picture where there 
is a random network of "fossil Stromgren trails" 
(Dupree & Raymond 1983) that ultimately might 
blend together. More extreme manifestations of 
such trails in denser media may have already been 
discovered by McCullough & Benjamin (2001) and 
Yagi et al. (2012), who observed faint, but straight 



and narrow lines of Ha emission in the sky (but 
were unable to identify the sources that created 
them). 

It is difficult to establish where in the contin- 
uum between the two extremes for the dispersal 
of ionization the true effects of embedded EUV 
and X-ray sources are to be found. For our treat- 
ment of the inffuence of these production sites 
for ionizing radiation, we will adopt the simplified 
premise that all of their photons are available to 
create a uniform but weak level of ionization every- 
where. This picture is not entirely correct, since 
one should expect that very near the sources some 
fraction of the ionizing photons are "wasted" by 
creating localized regions with much higher than 
usual levels of ionization. Such regions dissipate 
most of their ionization rapidly because their re- 
combination times are short. For this reason, our 
making use of a calculated average production rate 
of ionizing photons per unit volume will lead to an 
equilibrium equation that overestimates the space- 
and time-averaged level of ionization. (Later, it 
will be shown that this overestimate is of no real 
consequence because we obtain an answer that is 
still below that needed to explain the overall av- 
erage ionization level.) 

In the sections that follow, we consider three 
classes of sources that are randomly distributed 
throughout the neutral ISM and can in principle 
help to ionize it: main-sequence stars, active X- 
ray binaries, and WD stars. Luminous, early-type 
stars contribute large amounts of ionizing radia- 
tion, but most of their radiation completely ionizes 
the surrounding media and makes the gas virtu- 
ally invisible in the Ar I and O I lines. These stars 
also tend to be clustered inside the dense clouds 
of gas that led to their formation. 
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E (eV) 

Fig. 6. — Upper panel: External fluxes arising from the known X-ray background radiation (solid line) and 
a hypothetical, time-averaged flux (over a recombination time of « 1 Myr) from 3 supernova remnants at a 
distance of 100 pc (dashed line). The gray portion of the solid curve is in an energy range where the opacity 
to X-rays is high, and thus this radiation is not likely to penetrate into much of the gas that we observe. 
Lower panel: Internal average rates per unit volume for the injection of ionizing photons by the coronae 
of main-sequence stars (MS), active binaries (AB) and the photospheres of WD stars (WD). (A machine 
readable table of all 5 flux distributions vs. E shown in this figure is available in the online version of the 
Astrophysical Journal.) 
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The question may arise as to whether or not, 
by considering the embedded sources as a sepa- 
rate contribution that adds to the external radi- 
ation background discussed earher, we are possi- 
bly "double counting" some of the photons that 
could ionize the ISM. Kuntz & Snowden (2001) 
have computed the probable relative contribution 
of Galactic point sources that were not explicitly 
taken out of their measurements of the diffuse X- 
ray background, and they concluded that this con- 
tamination was, at most, only about 2 — 10% of the 
radiation that was thought not to arise from the 
Local Bubble. Within their highest energy bands 
(0.73 — 2.01 keV), they state that the contamina- 
tion could be as high as 51%. 

7.2.1. Main-sequence Stars 

For various kinds of point sources, the ratio of 
the emission of X-rays to photons in the V band 
is usually defined by the relation 

log(/x//y) = log fx + 0.4my + 5.37 , (17) 

where fx is the apparent X-ray flux of the source 
over a specified energy interval, expressed in 
erg cm~^s~-^, and my is its apparent visual mag- 
nitude0 The total output of X-rays from a source 
with an absolute magnitude My = should be 
5.11 X 10^''(/x//v) erg s~-^. Within each spectral 
class, there is a large dispersion in the measured 
values of log(/x//y), typically of order 1 dex, 
which is probably attributable to differences in 
stellar rotation velocities (Audard et al. 2000; 
Feigelson et al. 2004), variations in foreground 
absorption by the ISM, and time variability of the 
X-ray emission. While the most significant X-ray 
flares from stars can create spectacular increases 
in flux, their time-averaged effect has been esti- 
mated by Audard et al. (2000) to amount to only 
about 10% of the steady emission. On the basis of 
white-light monitoring of dwarf stars, Walkowicz 
et al. (2011) found that the duty cycle of flaring 
events is only of order a few percent. 

^Numerical values for this formula are not universal, since 
the adopted X-ray energy interval depends on which in- 
strument was used. For instance, surveys that used the 
Einstein Observatory , e.g. those reported by Maccacaro 
et al. (1988) or Stockc et al. (1991), used the pass- 
band 0.3 < E < 3.5 keV for fx, whereas those based 
on the ROSAT All-Sky Survey (RASS) (Krautter et al. 
1999; Agiieros et al. 2009) measured fx over the interval 
0.1 < B < 2.4 keV. 



For any given stellar spectral class with a char- 
acteristic absolute magnitude My that spans a 
range AMy along the main sequence and has 
luminosity function ^(My) stars mag~^pc~^, the 
energy density of X-rays per unit volume is given 

by 

F{E)ms = 1-74 X lO-2i0(My)AMy 
X ioi°g(/x//-)-o.4Afverg cm-^s-i . (18) 

If we use the mean values of \og{fx / fv) for differ- 
ent spectral classes listed by Agiieros et al. (2009) 
for the X-ray band 0.1 < i? < 2.4 keV, obtain 
values of My for these classes from Schmidt-Kaler 
(1982), define a luminosity function for stars in the 
disk of our Galaxy from the formula given by Bah- 
call & Soneira (1980), and then sum the results 
over all spectral types A-M, we obtain an average 
energy density equal to 1.90 x 10~^®erg cm~'^s~^. 
While we acknowledge that the spectral character 
of coronal emissions can vary for different stars 
along the main sequence, in the interest of sim- 
plicity we adopted for all cases a spectrum based 
on the differential emission measure (DEM) for the 
coronal emission from the quiet Sun, as defined in 
the CHIANTI database. To obtain a final photon 
emission rate per unit volume and energy in the 
ISM, we made use of this spectrum and normal- 
ized its energy output over the 0.1 — 2.4 keV band 
to the energy density factor given above. The re- 
sult is shown by the spectrum labeled "MS" in the 
lower panel of Figure [6l 

7.2.2. Active Binaries and Cataclysmic Vari- 
ables that Emit X-rays 

At high Galactic latitudes, most of the X-ray 
radiation above several keV originates from extra- 
galactic sources. Near the direction toward the 
Galactic center, however, Revnivtsev et al. (2009) 
found that at 4 keV about half of the X-ray emis- 
sion arises from point sources in the Galaxy that 
can be resolved by the Chandra X-ray Observa- 
tory, with the remainder coming from either a dif- 
fuse Galactic (hot gas) emission or an extragalac- 
tic contribution. To estimate the average emission 
per unit volume from these sources, we integrate 
the 2 — 10 kcV X-ray emission over the entire lu- 
minosity function 0(log L2_iokov)AB specified by 
Sazonov et al. (2006) to obtain the total energy 
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output 

L{2 - 10keV)AB 
/.34.0 

= / i2-lOkcV0(logi2-lOkcv)AB 
^22. 5 

X d(logL2-iokcv) = 2.13 X lO^Vg s-(l9) 

for all of the sources in our Galaxyll We may 
convert this value to an average volume emissiv- 
ity by multiplying it by the stellar mass density 
at our location (0.04 M0 pc~'^) divided by the 
total stellar mass of the Galaxy (7 x IO^^Mq), 
where both of these numbers were those adopted 
by Sazonov et al. (2006), and ultimately obtain 
a value 4.15 x 10~'^°erg cm^-^s^^. To define the 
spectral shape for the radiation emitted by these 
sources, we used the CHIANTI software to com- 
pute the emission from a plasma with an average 
of the DEM functions expressed by Sanz-Forcada 
et al. (2002, 2003jl that were constructed from 
their Extreme- Ultraviolet Explorer observations of 
various active binary sources. This spectrum was 
then normalized such that the flux in the 2—10 keV 
band matched the volume emissivity described 
above. The emission from active binaries that 
we derived F{E)ab is shown by the curve labeled 
"AB" in the lower panel of Figure [SI 

7.2.3. White Dwarf Stars 

WD stars that are hot enough to emit signif- 
icant fluxes in the EUV spectral range are much 
less numerous than the main-sequence stars con- 
sidered in Section 17.2.11 However their photo- 
spheres generate outputs in the EUV region that 
exceed by far the coronal emissions from indi- 
vidual main-sequence stars. This fact is demon- 
strated by the actual observations of the local 



'Note that the upper bound for the integration is at 
lO^^erg s~^. While the energy output over the whole 
Galaxy for neutron star and black hole X-ray binaries with 
individual outputs L2-iokcV > lO'^'^erg s~^ is substantial, 
these objects are so few in number that they can no longer 
be considered as embedded sources. The brightest such ob- 
ject in the sky, Sco X— 1, is at a distance of 2.8 kpc and cre- 
ates a flux at our location of only 2.8 X 10~^erg cm^'^s"^ 
over the 2 — lOkeV band (Grimm et al. 2002), which is 
substantially lower than the extragalactic background flux 
integrated over the whole sky. 

'These authors describe their DEM functions in terms of 
N{e)N{'a+)dV/d{logT), whereas the convention in CHI- 
ANTI assumes that the DEM is defined in terms of 
N(e)N{'R+)dV/dT. 



EUV sources compiled by Vallerga (1998), where 
he found that a vast majority of the detected 
objects were nearby hot WD stars0 Krzesin- 
ski et al. (2009) have measured the luminos- 
ity function for DA WDs 0(Mboi)wD, expressed 
in terms of (stars M^^jpc"^), in our part of the 
Galaxy using results from the Sloan Digital Sky 
Survey Data Release 4 (SDSS DR4) database. 
If we combine this information with the theo- 
retical computations of stellar atmosphere fluxes 
Fx (equal to 4 times the Eddington flux H\), 
expressed in the units erg cm~^s~^cm~^ (Ranch 
2003; Ranch et al. 2010), convert it to a phys- 
ical flux at the stellar surface J-{E) — 7.74 x 
10'^7rFxii;(eV)~^phot cm~^s~-^eV~\ and assume 
each star has a radius — O.O14r0 (Liebert et 
al. 1988), we can obtain a total emissivity per 
unit volume, 

FiE)wB = 47rr2(1.65 x lO^^^^) 

X J F{E)(t){M\,o\)wT>dMx,o\ phot cmT^s-^eYi^) 

A conversion from Mtoi to T^s is shown in the 
plot of the WD luminosity function presented by 
Krzesinski et al. (2009). The model fluxes for 
stars with Teff < 40, 000 K were assumed to arise 
from stars with pure hydrogen atmospheres, but 
stars with temperatures above this limit are known 
to have significant abundances of metals in their 
atmospheres because radiative levitation can over- 
come diffusive settling (Dupuis et al. 1995; Marsh 
et al. 1997; Schuh et al. 2002). Thus, for 
Teff = 50, 000 K and above, we used the fluxes 
for model atmospheres with [X] ~ \Y] = {) and 
[Z\ = — 1 , which have significant sources of opacity 
that reduce the radiation at energies E > 54 eV. 

As a check on the calculations described above, 
we can compute the X-ray energy outputs as a 
function of Toff over the 0.1 — 0.28 keV range, syn- 
thesize a luminosity function as a function of Lx 
in this band, and then compare the results to 
X-ray luminosity function derived by Fleming et 
al. (1996) from a RASS survey of WDs. Our 
source densities at the high end of the luminos- 
ity distribution compare favorably with the dis- 



The spectrum shown by Vallerga (1998) indicates that the 
my = 1.5 B2 II star Adhara (e CMa) dominates the local 
flux at energies below the He I ionization edge at 24.6 eV. 
This is an atypical situation, since N(H I) < lO^^cm"^ 
toward this star (Gry & Jenkins 2001). 
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tribution shown by Fleming et al., but we pre- 
dict a somewhat greater number of sources with 
Lx ^ lO^^erg s~^ due to a large space density of 
stars with 30, 000 < T^s < 40, 000 K. 

7.3. Interpretation of Internal Ionizations 

The treatments of the ionizations that arise 
from external and internal sources differ in a fun- 
damental way. External radiation above some en- 
ergy threshold is regarded as not being consumed 
by ionizing the gas, i.e., it is assumed to be unat- 
tenuated and thus each gas constituent is ionized 
at an appropriate rate F, as defined in Eq. [51 that 
is simply proportional to / a{E)F{E)dE. Here, 
<j{E) is the effective cross section for the combi- 
nation of the different ionization channels, as de- 
picted for Ar and H by the solid lines in Fig. [5l 

As indicated in the beginning of Section 17. 2[ 
for the internally generated radiation we switch to 
a very different concept and adopt the simplified 
premise that all of the photons emitted by embed- 
ded sources are used up by ionizing the various 
gas constituents that surround them. Thus the 
primary ionization rate for each kind of atom (or 
ion) X is given by 

Tp{X)^ j F{E)y{X,E)n{X)-^dE (21) 

where F{E) in this equation is the sum of all of 
the photon generation rates per unit volume spec- 
ified in Sections [7311 through [7X31 Each of the 
different species represented by X must compete 
with others for the photons that are consumed. 
Thus the equation includes a term y{X, E), which 
is a sharing function for the ionization rate that 
is represented by the relative probability that any 
photon with an energy E will interact with a given 
species X, 



yix,E) 



n{X)a{X,E) 
Yn(X')a{X',E) 

X' 



(22) 



where n{X) is the number density of X , <t{X, E) is 
the photoionization cross section of X at an energy 
E (likewise for X'), and the sum in the denomina- 
tor covers all of the major species competing for 
photons, i.e., X' = H", He°, He^, and heavy ele- 
ments whose inner shells respond to the more en- 
ergetic X-rays. For either the external or internal 
radiations, the secondary ionization rates V g and 



Fs' follow in proportion to the primary ones ac- 
cording to the descriptions given in Appendix \^ 
The ionizations from helium recombinations Fi 

±lc 

and Fj\g are treated as internal sources of ioniza- 
tion, and their rates are driven by the local densi- 
ties of He++, He+, and electrons, as described in 
Appendix [Bl 

7.4. Predicted Level of Ionization 

Given the computed rates of ionization in the 
previous sections, an evaluation of the electron 
fraction Xe = n(e)/n(Htot) will depend on both 
the temperature T and density n(Htot) for the 
gas. This fraction is higher than n(H+)/n(H°) 
because some electrons arise from the ionization 
of He. The coupling of the H and He ioniza- 
tion fractions is governed by Eqs. I13al through 
1151 Our observational constraint, which must ul- 
timately agree with the ionization calculations, 
arises from the values for [Ar I/O I], which respond 
to n(H+)/n(H'') in accord with Eq. [71 While the 
use of Pax in this equation will give an approxi- 
mate value for [Ar I/O I], a more accurate result 
emerges by replacing P^r by P^^, the derivation of 
which was described in Sections 16.41 and 16.51 Our 
goal will be to explore parameters that will match 
a computed value for [Ar I/O I] to the representa- 
tive observed value [Ar I/O I] -0.427 ± 0.11. 

We have no direct knowledge about the lo- 
cal values of n(Htot) that apply to the gas in 
front of the stars in this survey. We must there- 
fore rely on general estimates that have appeared 
in the literature. We can draw upon two re- 
sources. First, surveys of 21-cm emission indi- 
cate the amounts of H I in the disk at a Galacto- 
centric distance of the Sun, but difhculties in in- 
terpreting the outcomes arise from self-absorption 
effects and ambiguities in distinguishing between 
cold, dense clouds (CNM) and their surrounding 
WNM. Ferriere (2001) lists values for n(H°)wNM 
that are in the range 0.2 — 0.5 cm^'^. Dickey & 
Lockman (1990) state a value of 0.57 cm~'^, but 
it is not clear whether this excludes contribu- 
tions from the CNM. Kalberla & Kerp (2009) es- 
timate that the midplane density of the WNM is 
0.1 cm~'^, but this low value is difhcult to reconcile 
with measurements of the mean thermal pressures 
nT = 3800cm-3K by Jenkins & Tripp (2011) and 
a general recognition that Twnm < 10'* K. 
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Table 3 

Ionization Rates from Known Sources'* 



Ionization 


Rate 




Relative Contribution'^ 




Source'^ 


(10-"s-i) 




rs,H0 


rs.HcO 


rs,Ho+ 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


r(H'')oxt. 


5.26 


0.20 


0.22 


0.52 


0.013 


0.048 


r(H°)Ms 


6.72 


0.26 


0.044 


0.57 


0.11 


0.012 


r(HO)AB 


0.649 


0.035 


0.015 


0.52 


0.22 


0.21 


r(H")wD 


2.19 


0.92 


0.030 


0.044 


6.3e-3 


9.4e-6 


rffe+(H") 


0.166 












rHeo(HO) 


4.94 












Ccr(HO) 


12.3 












Total r(H°) 


32.3 












r(He°)ext. 


31.4 


0.87 


0.035 


0.084 


2.2e-3 


7.5e-3 


r(He°)MS 


25.7 


0.80 


0.012 


0.15 


0.028 


3.0e-3 


r(He°)AB 


0.861 


0.33 


0.011 


0.36 


0.16 


0.14 


r(He°)wD 


17.0 


0.99 


4.5e-3 


5.9e-3 


8.6e-4 


1.2e-6 


rHe+(He°) 


0.916 












rHeo(Hc°) 


7.10 












CcR(He°) 


14.4 












Total r(He°) 


97.3 












r(He+)ext. 


16.9 












r(Hc+)Ms 


8.32 












r(He+)AB 


0.0978 












r(He+)wD 


0.224 












THe+ (HC+) 


0.317 












Total r(lle+)'^ 


25.9 












r(Ar°)ext. 


251. 


0.94 


0.018 


0.042 


l.le-3 


3.9e-3 


r(Ar°)MS 


82.2 


0.77 


0.014 


0.18 


0.035 


3.9e-3 


r(ArO)AB 


10.5 


0.77 


3.6C-3 


0.12 


0.053 


0.049 


r(ArO)wD 


45.3 


0.99 


5.6e-3 


8.2e-3 


1.2e-3 


1.8e-6 


rffe+(ArO) 


1.35 












r^,o(Ar°) 


93.6 












CcR(Ar°) 


113. 












Total r(Ar°) 


597. 












r(Ar+)ext. 


123. 












r(Ar+)MS 


47.2 












r(Ar+)AB 


1.57 












r(Ar+)wD 


49.0 












r//e+(Ar+) 


0.969 












Total r(Ar+)d 


222. 
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'^The internal rates (with subscripts MS, AB and WD) are based on a 
volume density n(Htot) = 0.5 cm~"^. Such rates scale inversely with den- 
sity (although not exactly so, because the secondary ionization efficiencies 
change when Xe changes). 

^Fractions of the values shown in column (2). See Egs. [51 \Xl\ IA21 [X3l 
and Ell 

'^Meaning of the subscripts that follow the different forms of F: 
ext. — external X-ray background radiation (Section 17. MS = embed- 
ded main-sequence stars (Section I7.2.ip . AB = embedded active binary 
stars (Section [7X2]) ■ and WD = embedded WD stars (Section 17X31) . 

'^The total ionization rates for the ions He+ and Ar+ do not include the 
ionizations from secondary electrons or cosmic rays. Hence these totals 
underestimate the true rates. 
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Heiles & Troland (2003) estimate that an over- 
all average (n(H")wNM)) = 0.28 cm~'^ translates 
into local densities of 0.56 cm~'^ if the WNM has 
a volume filling factor of 0.5. A second method 
for estimating n(H°)wNM is to note where the 
WNM branch of the theoretical thermal equilib- 
rium curve of Wolfire et al. (2003) intersects the 
average thermal pressure p/k — 3800cm~''K of 
Jenkins & Tripp (2011). This exercise yields a 
value of 0.4 cm~'^. 

Table [3] shows a detailed accounting of the ion- 
ization rates from the different radiation sources, 
under the condition that the value of Xe corre- 
sponds to what the equilibrium equations yield 
for a WNM at T = 7000 K. The entries in this 
table reveal that for neutral hydrogen the ioniza- 
tion caused by all of the sources of EUV and X-ray 
photons result in a combined total rate r(H'') that 
is about 1.6 times the rate of ionization by cosmic 
rays CcR- For neutral He and Ar, the photon ion- 
ization rates are several times higher than those 
from cosmic rays. The changing patterns in the 
distribution of different ionization mechanisms re- 
vealed in Columns (3) to (7) in the table reflect 
differences in the distribution of fluxes with en- 
ergy shown in Fig. [B) For sources with hard spec- 
tra, about half of the hydrogen ionization arises 
from the ionization of He°, which produces en- 
ergetic electrons that can cause secondary coUi- 
sional ionizations, i.e., the process associated with 

fj(,o (H"). It is only for the very soft spectrum 
of the collective radiation from WD stars that we 
find that Fp(H°) strongly dominates over the other 
ionization routes. 

The ionization rates and concentrations of vari- 
ous primary constituents are coupled to each other 
by the network of reactions described in Section |6] 
and Appendices |^ through The first column 
of Table m lists various quantities of interest, and 
their derived values under the assumption that 
n(Htot) = 0.50 cm-3 and T = 7000 K are given 
in Column (2). (The remaining two columns of 
this table will be discussed later in Section |8l) 

Figure [7] shows the predicted outcomes for 
[Ar I/O I] after we perform the calculations, again 
using the equations and reaction rates given in 
Section [S] and the Ar ionization rates derived 
from the information in Sections 17.11 and 17.21 
The upper curve in this figure that represents 
'^(Htot) — 0.50 cm^"^ is clearly inconsistent with 



our findings for [Ar I/O I]. At T = 7000 K, a lower 
value ?i(IItot) — 0.14cm~'^ is consistent with the 
upper error bound for [Ar I/O I]. At lower temper- 
atures, however, the predictions for [Ar I/O I] once 
again are found to be considerably above the ob- 
served values. It is not until a density of 0.09 cm""^ 
is reached, a value that is unacceptably low, that 
the predicted ionization conditions for T = 7000 K 
fit comfortably with the nominal value from the 
observations. Even here, however, this result is 
not fully satisfactory because a good fraction of 
the WNM is known to be at temperatures well 
below 7000 K, where the gas is thermally unstable 
(Heiles & Troland 2003). 

Our overall conclusion is that the ionization 
rates arising from what we consider to be known 
sources of ionizing radiation are not able to main- 
tain a level of ionization in the WNM that is con- 
sistent with our low observed values for [Ar I/O I]. 
Our quest to resolve this problem by exploring 
some possible supplemental means for ionizing the 
medium will be addressed later in Section [8l 

1.0 I I I I ITT1 1 1 1 — I I I I I I 1 1 1 — I I I I I 
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100 1000 10000 

T(K) 

Fig. 7. — Calculated values for lol^''^/^^! for three 
different densities n(Htot) and a range of differ- 
ent gas temperatures T shown on the x axis, com- 
puted using estimates for the ionization rates from 
known sources of radiation, both external and in- 
ternal. A general value for lOl^'''^/*-''! that arose 
from the survey is indicated by the horizontal line, 
and the shaded band shows the range of possible 
systematic errors. 
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Table 4 
Gas Properties'" 





Without SNR 


With SNR 


Greater Low 


Quantity 


Contrib.'" 


Contrib.'^ 


Energy Penetration*^ 


(1) 


(2) 


(3) 


(4) 


PAr" 


13.1 


19.3 


12.6 


P' f 


13.4 


22.8 


15.9 


n(e) (cni~^) 


0.018 


0.033 


0.052 


n(HO) (cm-3) 


0.48 


0.48 


0.46 


n(H+) (cm-3) 


0.015 


0.025 


0.036 


n(He°) (cni-3) 


0.047 


0.042 


0.035 


n(He"'") (cm^^) 


3.0e-03 


8.1e-03 


0.014 


n(He "'""'") (cm~^) 


2.1e-05 


2.1e-04 


l.Oe-03 


Total r(HO) (s-i) 


3.2e-16 


9.2e-16 


2.1e-15 


4 (10-"erg s-^H atom-^)s 


4.3 


6.2 


8.5 



'"For n(Htot) = 0.50 cm-^ and T = 7000 K. 

'^Apphes to the photoionization rates expressed in Table [3] from the known sources 
only. 

'^The radiation from known sources is supplemented by the average flux over 1 Myr 
from 3 supernova remnants created by explosions with energies E'en = 3 x 10^" erg in 
a medium with n(H) = 1.0 cm~^ located at a distance of 100 pc. This supplemental 
radiation is shown by the dashed line in the upper panel of Fig. |6l 

Allows for a low energy cutoff of 60 eV instead of 90 eV for the external radiation 
background; see Section [5^ 

''As defined in Eq. [5] 

*As defined in Eq. \W[ 

^Cooling rate from radiation by C+ in its excited fine-structure state. 
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7.5. Dependence of Recombination with 
Time 

An argument that helps to support the concept 
that the internal sources spread their ionization 
rather evenly throughout the medium is that they 
move at a rate that makes their local residence 
time short compared to a characteristic e-folding 
time <iccomb. for the decay of the proton density 
from an initial high value to some end equilibrium 
state n(H+)eq. as t -> oo, 



(-rccomb. 



n(H+) - n(H+) 



cq. 



dn{}i+)/dt 



(23) 



where 



dnOi+)/dt^~ a(H°,T)n(e) 
+ ag(H",n(e),G,T)n(Htot) 
+ r(H")ln(H+) + r(H")n(Htot) 



(24) 



This expression overlooks the complications aris- 
ing from charge exchange reactions, and when we 
evaluate the trend of n(H"'") with time we assume 
that n(e) is always equal to 1.2n(H+), as is the 
case when the gas is in a steady-state ionization 
condition. Another simplification is that T, and 
hence a(H°, T), remains constantly When the gas 
is highly ionized, r(H°) is not equal to the value 
that we computed for n(H+)oq. — 0.015 cm~^ be- 
cause the secondary ionization processes change 
with n(e) and depend on the strength of the he- 
lium ionization. However, we can ignore this com- 
plication that occurs at early times because the 
recombination terms in Eq. [24] are considerably 
larger than the terms involving r(H'^). The ion- 
ization rate influences the character of the de- 
cay only when the gas is weakly ionized. Hence 
there is no harm in declaring that at all times 
r(H°) = 3.23 X lO-i^s-i, as stated in Table El 
We can adopt for the final state (equilibrium) 
densities the values ?i(H+)cq. = 0.015 cm~'^ and 
n(e)eq. = 0.018 cm-3 i^^^^^i Column (2) of Ta- 
ble H 



^^A lack of variation in T is probably a safe assumption for 
isobaric recombination, but in the isochoric case T may de- 
viate to lower values at intermediate times (Dong & Draine 
2011), and this would increase the value of a(H°,T) and 
make the recombination more rapid. This could be of im- 
portance for very large regions of space ionized by the ra- 
diation from old SNRs, which will be considered later in 
Section [Q 



The change of irocomb. with time and the relax- 
ation of n(H+) from a fully ionized condition to 
n(H+)oq. for n(Htot) = 0.50 cm'^ and T = 7000 K 
is illustrated in Fig. [51 Values for trocomb. start at 
0.13 Myr for an initial ri,(H+) = 0.50 cm~'^ and in- 
crease to 1.0 Myr when n(H+) reaches 0.036 cm~^ 
ai t — 1.77 Myr, beyond which irocomb. very grad- 
ually climbs toward a steady value of 1.68 Myr and 
the subsequent decay in n(H+) toward n(H+)oq. = 
0.015 cm~^ is almost purely exponential. 



(0 

> 




2 4 

f[Myr] 



Fig. 8. — The behavior of irccomb. (upper solid 
curve) and n(H+) (lower solid curve) with time t, 
starting from a fully ionized condition to an equi- 
librium state ri(H-|-)oq. for n(Htot) — 0.50 cm~^ 
and T = 7000 K. Values of ^rocomb. are defined by 
Eq. [23l and n(H-l-) as a function of t can be found 
by solving the differential equation for dn{ll+) / dt 
given in Eq. [24l By solving a similar equation for 
Ar"*" and comparing it to that for H"*", we obtain 
the trend for [Ar I/H I] shown by the dashed line. 

Of the three different classes of internal sources, 
WD stars emit the softest radiation, which means 
the infiuences of their ionizations can be more 
sharply bounded than those of the others. While 
this consideration may present a challenge to our 
assumption about the uniformity of ionization in 
space, we note that these stars have observed ra- 
dial velocity dispersions of about 25 km s~^ (Fal- 
con et al. 2010), or rms speeds of 43 km s~^, which 
are consistent with the transverse velocity disper- 
sion measurements reported by Wegg & Phinney 
(2012). This indicates that WDs typically move 
about 76 pc during the time interval that n(H+) 
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decays from 0.50 to 0.036 cm~'^. This amount of 
travel is considerably larger than the radius of a 
Stromgren Sphere 

" n(H)2a(H0,r) ) ^^^> 

that would be established if a WD star were sta- 
tionary. (Recall that F(i?)wD is the flux emitted 
at the surface of the star with a radius .) A typi- 
cal value for rs in the WNM would be about 1.8 pc; 
this value applies to a WD star with r, — O.O14r0 
that has Tcff — 50, 000 K and is situated within a 
gas with a density n(Htot) — 0.50 cm^'^. 

One can view the dispersal of ionization caused 
by the star's motion in the context of the sim- 
plified description of "cometary H II regions" 
described by Raga (1986); the dimensionless 
elongation parameter in this development /3 = 
3v/[n{litot)a(lf,T)rs] = 14 if rs = 1.8 pc and 
V = 43 km s~^, which would result in a tube of 
ionization with a radius ~ 0.5rs and a very long 
tail. The fact that the ionizations produced by 
WD stars could be highly diluted because they 
are distributed over large volumes could explain 
why a sensitive survey of Ha emission from regions 
around such stars carried out by Reynolds (1987) 
yielded only one detection out of nine targeted 
regions surrounding these stars. 

8. PROPOSALS FOR ADDITIONAL PHO- 
TOIONIZATION 

Given that the outcome for a reasonable value 
of n(Htot) and the ionization rate from known 
sources discussed in Section [7] do not agree with 
our observations, we must explore some alterna- 
tives for boosting the average rate of photoioniza- 
tion in the WNM. Two possibilities are discussed 
in the following sections. 

8.1. Ionization from Previous Supernova 
Remnants 

A popular theme in the ISM literature of early 
1970's was the consideration that the low density 
portions of the ISM were ionized impulsively by 
UV and X-ray radiation from supernovae and their 
remnants (SNRs), which represented sources of 
ionization that had limited durations and that ma- 
terialized at random locations and times (Bottcher 



et al. 1970; Werner et al. 1970; Jura & Dalgarno 
1972; Gerola et al. 1973; Schwarz 1973). A par- 
ticularly instructive example that shows the wide 
variations in temperature and fractional ionization 
can be seen in the results from a numerical sim- 
ulation by Gerola et al. (1974). This topic has 
also been approached by Slavin et al. (2000), who 
estimated the effects of irradiation of the neutral 
medium by SNRs, with a special emphasis on its 
influence on gas in the Galactic halo. The spec- 
tral character and strength of the emitted radia- 
tion from a remnant depends on a combination of 
supernova energy and the density of the ambient 
medium (Mansfield & Salpeter 1974). The emis- 
sion of soft X-rays as a remnant evolves can last 
up to a late phase when radiative cooling causes a 
dense shell to form, which could be opaque to the 
lowest energy X-rays. However some of the radia- 
tion may continue to escape afterward if instabil- 
ities cause the shell to break up and create gaps 
(Vishniac 1994; Blondin et al. 1998). The magni- 
tudes of such instabilities are uncertain, since they 
might be diminished if either the ambient medium 
has a low density and high temperature (Mac Low 
& Norman 1993) or if the shock is partly stabilized 
by an embedded transverse magnetic field (Toth 
& Draine 1993). 

Chevalier (1974) has computed the strength 
and distribution over energy for the radiation 
emitted by various models for SNRs. As a rep- 
resentative example, we can use the emission over 
a time interval of 2.5 x 10^ yr from his model that 
had an energy Esn = 3 x 10^" erg and was ex- 
panding into a medium with an average density 
of Icm""^. If we assume that there was a suc- 
cession of about 3 such remnants that occurred 
once every Myr (i.e., of order trccomb. defined in 
Eq. [23]for n(H+) = 0.036 cm-^; see Fig. [8]) and 
they were at a distance of 100 pc from some loca- 
tion in the ISM, the resulting time-averaged flux 
shown by the dashed line in the upper panel of 
Fig.[n]would increase the total ionization of H from 
r(H°) = 3.2 X lO^i^s^i for the steady ionization 
sources to the much higher value of 9.2 x 10~^^s~^ 
and create conditions that would yield values for 
[Ar I/O I] that are virtually the same as those de- 
picted by the curve labeled n(Htot) = 0.14 cm~^ 
in Fig. [71 but in this case for a density n(Htot) = 
0.50 cm-3 At T = 7000 K, the densities of var- 
ious constituents are represented by the numbers 
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that are shown in Column (3) of Table 21 For this 
temperature, the calculation of [Ar I/O I] is just 
consistent with the upper error bound from the 
observations. 

A stronger flux from SNRs would be needed to 
make the calculated value for [Ar I/O I] match 
nominal value for the observed ones. To obtain 
a curve that matches the lowest curve (labeled 
'^(Htot) — 0.09 cm^'^) in Fig. [7] and still maintain 
an assumed local density 7i(Htot) — 0.50 cm^'^, 
we would require either an increase in the aver- 
age energy of the supernovae (SNe), or the rate 
would need to be increased from 3 to 8 SN Myr^^ 
at a distance of 100 pc. In this circumstance, 
n(e) = 0.047 cm"3. 

While one might be tempted to ask whether 
some average SN rate in our location within the 
Galaxy (van den Bergh & McClure 1994; McKee & 
Williams 1997; Ferriere 2001) is consistent with at 
least 3 SNe Myr~^ at separations of about 100 pc 
from some representative location, it is probably 
more realistic to draw upon actual estimates of the 
recent history of explosions in the general neigh- 
borhood of the Sun, which have been inspired by 
the evidence of how the ISM has been disrupted 
by the expansions of the explosion remnants. The 
Local Bubble (see footnote 2 in Section [2]) and the 
neighboring radio emitting Loop I superbubbld^. 
which intersect one another (Egger & Aschenbach 
1995), are both thought to have been created by 
a series of SN explosions that occurred over the 
past 10-15 Myr (Maiz-Apellaniz 2001; Berghofer 
& Breitschwerdt 2002). Fuchs et al. (2006) an- 
alyzed the numbers and mass functions of stars 
within nearby 0-B associations and traced their 
motions back in time. From this information, they 
concluded that the Local Bubble was created by 
14-20 SNe. This estimate is consistent with the 19 
SNe that Breitschwerdt & de Avillez (2006) used 
in their hydrodynamic simulation of the creation 
of the Local Bubble. 

Within a shorter time frame, there is some evi- 
dence that atoms from a nearby SN were deposited 
on the Earth some 2.8 Myr ago, as revealed by an 
enhancement of ^^Fe in a thin layer within deep- 
sea ferromanganese crust (Knie et al. 2004). Fi- 



Loop I is coincident with tiie Nortii Polar Spur of X-ray 
radiation. Its center is estimated to be only 180 pc away 
from us (Bingham 1967). 



toussi et al. (2008) attempted to replicate this 
result in a marine sediment, but the outcome did 
not agree with the earlier result at the same level of 
significance. However, there are alternative means 
for investigating past depositions of SN elements 
in terrestrial samples (Bishop & Egli 2011; Feige 
et al. 2012), and new results may emerge soon. 

If we add to the Local Bubble contribution the 
past SNe in the Sco-Cen association that created 
Loop I (Iwan 1980; Egger 1998; Diehl et al. 2010), 
one can imagine that, to within the uncertainties 
in the X-ray production rates, the time-averaged 
level of ionizing radiation could conceivably be of 
the same order of magnitude as that discussed in 
the preceding paragraphs. 

A legitimate question to pose is whether or not, 
at some intermediate stage of recombination, the 
value of [Ar I/HI] makes a brief excursion to a 
level below the final equilibrium state. If this 
were the case, we might be able to relax the re- 
quirement for a high frequency of SN explosions. 
The recombination rates Q;(H°,r) and a(Ar'') are 
nearly equal to each other, but the dust grain re- 
combination rate for Ar should be much less 
than that for H because the thermal velocities 
of Ar are lower. Thus, at early times where re- 
combinations still dominate over the steady-state 
ionization rate but the electron density has de- 
creased to a point that neutralization of Ar and H 
by collisions with dust grains become important, 
one could imagine that the difference in grain re- 
combination rates might be influential in allowing 
the H+ to recombine more quickly than Ar"*^. In 
order to investigate this possibility, we can evalu- 
ate the time history of the Ar recombination using 
Eqs. [231 and [Ml but with a substitution of "Ar" for 
each appearance of "H." A resulting comparison 
of the Ar recombinations to those of H shows that 
[Ar I/H I] shows a steady decrease from an ini- 
tial value slightly greater than 1.0 (when the gas 
is just starting to recombine) to the final equilib- 
rium state for this quantity without any excursion 
to lower values. This behavior is illustrated by the 
dashed line in Fig. [H 

8.2. Penetration of the WNM by Low- 
energy Photons 

Until now, we have regarded the emission of X- 
rays from the hot (T ~ 10^ K) gas in our Galaxy 
as an external source of ionizing radiation that 
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must penetrate the bulk of the WNM regions un- 
der study. With this condition in mind, we es- 
tabhshed a 90 eV cutoff for the external radia- 
tion, below which the X-rays were assumed to be 
absorbed in the outer layers of the neutral gas. 
However, some of the soft X-ray emitting gas may 
be threaded through the WNM, a picture that is 
reminiscent of the network of hot interstellar tun- 
nels proposed by Cox & Smith (1974). Hot, inter- 
spersed gas could allow some or all of the neutral 
medium to have access to lower energy radiation. 
Observations indicate that O VI, an indicator of 
gas that is coUisionally ionized at T ^ 3 x 10^ K, 
can be seen in emission (Dixon et al. 2006; Otte 
& Dixon 2006) and absorption (Jenkins 1978a, b; 
Bowen et al. 2008) almost everywhere in the disk 
of the Galaxy. The simulations of SNe explod- 
ing at random in the Galactic disk by Ferriere 
(1995) and de Avillez & Breitschwerdt (2005a, b) 
indicate that this hot gas should be sufficiently 
pervasive and frothy that, even though our mea- 
surements penetrate through column densities well 
in excess of lO^^cm"^, individual parcels of gas 
may be exposed to adjacent sources of radiation 
through significantly lower column densities. In- 
deed, there could be unusually strong X-ray emis- 
sion at the boundaries where the neutral regions 
come into contact with the hot gas, where charge 
exchange reactions between neutrals and highly 
ionized species could produce an enhancement of 
soft X-ray emission over that which is emitted by 
just the hot gas (Wang & Liu 2012). 

It is clear from the gray parts of the exter- 
nal radiation spectrum depicted in Fig. [6] that 
even a modest lowering of the cutoff below our 
adopted 90 eV threshold will create a significant 
increase in the ionization of the WNM. For in- 
stance, if the threshold were dropped to 60 eV, an 
energy just below the strong M-line complex of 
emission features from highly ionized Fe at 70 eV, 
the level of ionization would increase to an ex- 
traordinarily large total hydrogen ionization rate 
r(H°) = 2.1 X lO^i^s-i. For n(Htot) = 0.50 cm'^ 
and T = 7000 K, this increase would create a value 
for [Ar I/O I] that is consistent with the upper 
error bound for the measured outcomes. In this 
particular case, most of the ionization is caused 
by photons with energies near 70 eV. For this en- 
ergy, Pai{E) is near its minimum value (see the 
lower panel of Fig. [5] and the entry in Column (4) 



of Table H]). As a consequence, in order to match 
our observations, Eq. [7] shows us that the required 
n(H+) (and thus r(H°)) must be increased beyond 
the value needed from the time-averaged SNR ra- 
diation with its larger value of Pat- 

A consideration that disfavors the presence of 
a strong ionizing flux at energies in the vicinity of 
70 eV, at least at our location, is that it does not 
seem to be present at anywhere near the intensity 
level shown in Fig. [51 For instance, from observa- 
tions by the Cosmic Hot Interstellar Plasma Spec- 
trometer (CHIPS) Hurwitz et al. (2005) stated an 
upper limit for the flux emitted by the Fe lines near 
70 eV. This limit was consistent with an emission 
measure for a plasma with a solar abundance pat- 
tern at T = 10^ K that is less than 5% of the value 
EM = 10^^ '^^cm^^ that we used for reconstruct- 
ing the emission spectrum from the soft compo- 
nent of the X-ray background described in Sec- 
tion [73] Jelinsky et al. (1995) and Bloch et al. 
(2002) similarly found flux upper limits at low en- 
ergies (but not as stringent as those from CHIPS) 
that were much lower than model predictions at 
this energy that we obtained from fits to the ob- 
served radiation at higher energies. A marginal 
detection of the Fe emission complex toward high 
Galactic latitudes by McCammon et al. (2002), 
FiE) = 100 ± 50 phot cm^^s^^sr-i is also well 
below the peak at 70 eV that appears in our recon- 
structed flux. The weakness of these fluxes could 
be caused either by a deficiency of Fe below the 
solar abundance ratio in the hot emitting gases or 
by our inability to see hot gas regions whose low 
energy radiation is not absorbed by intervening 
neutral material. 

9. DISCUSSION 

We have now reached a point where it is appro- 
priate to investigate the consistency of the newly 
derived WNM ionization levels with other obser- 
vational and theoretical findings, along with some 
consequences of our results on various relevant 
physical processes. 

9.1. Other UV Observations 

Along a number of sight lines, there have been 
detailed investigations of UV absorption lines ob- 
served at high enough velocity resolutions to iden- 
tify which components came from either mostly 
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neutral or mostly ionized clouds (Spitzer & Fitz- 
patrick 1993, 1995; Fitzpatrick & Spitzer 1994, 
1997; Wood & Linsky 1997; Holberg et al. 1999; 
Welty et al. 1999; Jenkins et al. 2000a; Gry & 
Jenkins 2001; Sonnentrucker et al. 2002). These 
investigations relied either on the relative popula- 
tions of the excited fine-structure level of C+ or 
the ratios of ions to neutrals for various elements. 
The ion-to-neutral ratios of different elements gave 
very different outcomes for n{e) (but systemati- 
cally went up and down together from one velocity 
component to the next); these disparities probably 
arose as a result of a lack of a good understand- 
ing of the physical processes involved. Welty et 
al. (2003) have presented a good summary of the 
results that exhibited this problem. 

Values of n{e) identified with the WNM by the 
investigators cited above generally ranged between 
0.04 and 0.12 cm^^, which seem to be higher than 
our values given in Table SI A common theme 
in the discussions of these results was that such 
high partial ionizations could not be attributed to 
the known EUV, X-ray and cosmic ray ionization 
rates - a conclusion that is stated once again from 
our [Ar I/O I] results. 

Electron densities for clouds embedded in the 
Local Bubble have an average value of 0.11 cm~^ 
(Redfield & Falcon 2008), a result that once again 
relied on comparisons of excited C+ to other 
species (either unexcited 0+ or S+). The frac- 
tional ionization of these clouds appears to be 
much higher than what we found for the general 
WNM outside the Local Bubble, if one assumes 
that the characteristic total density ri(Htot) ~ 
0.2 cm~'^ within the clouds (Redfield & Linsky 
2008). By solving for the time-dependent ioniza- 
tion in Eq. [24] for n(Htot) — 0.2 cm^'^, we find 
that the decay from a fully ionized condition to 
an approximately half ionized state takes only 
0.40 Myr. It is likely that this higher level of ion- 
ization might be explained by supplemental radi- 
ation from evaporative boundaries that surround 
the clouds (Slavin & Frisch 2002), if indeed they 
are embedded in a very hot, low density gas, or 
by the infiltration of some ionizing radiation from 
e CMa, which strongly dominates over other radi- 
ation sources within the Local Bubble at low en- 
ergies (Vallerga & Welsh 1995). 



9.2. Pulsar Dispersion Measures 

Cordes & Lazio (2002) have derived character- 
istic electron densities for three volumes within 
about 1 kpc of the Sun. Inside the Local Bub- 
ble, their model seems to best fit an average den- 
sity of only {n{e)) = 0.005 cm~^, which is not sur- 
prising because much of the volume has probably 
been cleared of material by SN explosions, except 
for some isolated, warm clouds with an average 
filling factor in the range 5.5 — 19% (Redfield & 
Linsky 2008). Regions outside the Local Bubble 
have higher densities: two volumes that they stud- 
ied yield (n(e)) = 0.012 and 0.016 cm-3_ how- 
ever, these regions are identified with ellipsoidal 
volumes designated as either a "local superbub- 
ble" (LSB) or a "low density region" (LDR), so 
they too could have significant voids that would 
dilute the apparent electron densities. The values 
of (f^(e)) given for these regions could be consis- 
tent with our measures of n(e) given in the last two 
columns of Table |4] for the WNM if this medium 
had a filling factor of about 1/3 in these regions 
and there were no significant contributions from 
the warm ionized medium (WIM). 

An independent analysis of dispersion measures 
was carried out by de Avillez et al. (2012) for 
24 pulsars with known distances between 0.2 and 
8 kpc from the Sun and with \z\ < 0.2 kpc. They 
found a distribution of n(e) outcomes that was 
consistent with a log-normal distribution centered 
on log ri(e) = —1.47 and a dispersion cr[logri(e)] = 
0.17. The fact that their representative values of 
n(e) are higher than the determination of Cordes 
& Lazio (2002) and close to our findings based on 
[Ar I/O I] may be accidental, since their results 
may be strongly influenced by sight-line intercep- 
tions of fully ionized regions. 

9.3. Emission Lines 

Observations of diffuse line emissions in the sky 
most generally apply to probing the physical con- 
ditions in fully ionized gases, either the bright H II 
regions around hot stars or the WIM (Reynolds et 
al. 1977, 2002; Haffner et al. 1999; Madsen et al. 
2006). Most of the emissions are dominated by 
contributions from species that are expected to be 
abundant in such highly ionized media, and they 
should overwhelm any contributions from the very 
dilute ionization in the WNM. Nevertheless, there 
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are a few cases where emission line fluxes from 
some neutral species in the WIM are expected to 
be very weak, and thus a contribution from the 
WNM might in principle be identified. We need 
to explore whether or not the predictions for line 
strengths from a medium with our enhanced elec- 
tron densities are not in serious violation of detec- 
tions or the upper limits for the fluxes. We explore 
three such cases in the subsections that follow. 
The first test involving line emission from He re- 
combinations is especially important, because we 
predict that the fraction of singly-ionized He could 
be as high as 16—28% of the total amount of he- 
lium. 

9.3.1. He \5876 Recombination Radiation 

Reynolds & Tufte (1995) attempted to compare 
the strength of the recombination line of He° at 
5876 A to that of Ha in parts of the sky away 
from well defined H II regions. Their motive was 
to determine the hardness of the radiation that 
maintains the ionization in the WIM. While an 
explicit upper limit for the He recombination line 
flux was not stated by Reynolds & Tufte (1995), 
we estimate that the flux in the two directions 
that they sampled was found to be less than 0.1 R 
(R = Rayleigh = 106/(47r)phot cm"^sr-is-i). 

Using the line emission rates given by Benjamin 
et al. (1999) we find that for n(e) = 0.04 cm-^, 
n(He+) = 0.01 cm-3, and T = 7000 K (represen- 
tative values for the enhanced ionization cases pre- 
sented in the last two columns of Table the 
emission should be 3.0 x lO^^^phot cm s . If 
we were to propose that the emission is seen over 
a path of 200 pc with no extinction (i.e., imagine a 
length of 400 pc with a filling factor for the WNM 
of 50%), we could expect to find an emission of 
1500 phot cm"^sr"is-i = 0.02R. This expecta- 
tion is well below the sensitivity of the observa- 
tions by Reynolds & Tufte (1995), so our predicted 
intensity from the enhanced electron density and 
ionization of He does not violate their upper limit. 

9.3.2. Emission of [0 I] \6300 from Electron 
Collisions 

In order to determine the neutral fraction of 
H in the WIM, Reynolds et al. (1998) measured 
the strength of the [O I] A6300 line emitted in 
parts of the sky that had a uniform, moderately 



strong Ha emission (but away from obvious H II 
regions excited by stars), much as they had done 
for the He recombination line discussed above. In 
three different directions, they detected intensities 
of 0.2, 0.09, and 0.11 R. 

According to Federman & Shipsey (1983), elec- 
tron collisions dominate over those by hydrogen 
for the excitation of the ^I?2 state of neutral oxy- 
gen when the electron fraction exceeds 1.5 x 10^^. 
Hence, we can ignore hydrogen impacts. From the 
fitting formula of Pequignot (1990) to the colli- 
sion strengths computed by Berrington & Burke 
(1981), we derive a collisional rate constant Ce for 
electron excitations to be 8.0 x 10~^^cm'^s~^ at 
T = 7000 K. If we again make the conservative as- 
sumption expressed in Section r4.2.1l that at the low 
densities of the WNM the depletion of O is negli- 
gible, we expect that n{0°) = 2.7 x 10~^cm~^ if 
n(H°) = 0.47 cm"^. For n{e) — 0.04 cm^^, we ex- 
pect the emissivity to be equal to n{0'^)n{e)Ce 
multiplied by a branching fraction 0.76 (Froese 
Fischer & Tachiev 2010) for the proportion of de- 
cays from the level to the lower level. This 
product equals 6.6 x 10~^®phot cm~'^s~"'^, which 
should produce 0.41 R over a path of 200 pc. This 
value is greater than the three measurements by 
Reynolds et al. (1998). The magnitude of this vi- 
olation is not large, considering the uncertainties 
of our assumptions about the lack of depletion of 
O and the adopted path length estimate. Also, 
for temperatures less than 7000 K, the expected 
strength of the emission will be considerably less: 
for instance, at T = 5000 K, the emission should 
be 4.5 times weaker than at 7000 K. 

9.3.3. Emission of [N I] X5201 from Electron 
Collisions 

An upper limit of 0.13 R for the [N I] A5201 
line was determined for a single direction in the 
sky by Reynolds et al. (1977). If we adopt the 
same calculations as in the previous section for 
O I, but make the substitution that n(N°) = 8.0 x 
10~^n(H°) = 3.8 X 10~^cm~^ (again, assuming 
no depletion), and Ce — 8.6 x lO^^^cm'^s"^ at 
T = 7000 K (Tayal 2006), we obtain an emissivity 
equal to 1.3 x 10~^^phot cm^'^s"^, which yields an 
intensity of 0.08 R over a 200 pc path. This value 
is below the upper limit determined by Reynolds 
et al. (1977). 
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9.4. Cooling Rates from Carbon Ions and 
Oxygen Atoms 

A major coolant for the neutral ISM is the 
singly charged carbon ion (Dalgarno & McCray 
1972; Wolfire et al. 1995, 2003), whose ex- 
cited fine-structure level in the ground electronic 
state can be excited by collisions with electrons 
and hydrogen atoms. The cross section for exci- 
tation by electrons is substantially greater than 
that for neutral hydrogen. For this reason, the 
level of partial ionization is an important factor 
in the excitation rate. After excitation, an energy 
loss can occur because the excited level can un- 
dergo a spontaneous radiative decay with a tran- 
sition probability A21 = 2.29 x 10~^s~^ to the 
lower state (Nussbaumer & Storey 1981), a 

process that liberates a photon with a wavelength 
of 158 /im. If we assume that the abundance ratio 
(C+/HO) = 9.5x10-5 (see footnote 1 inSectionP, 
we can use the densities that we derived for n(H°) 
and n(e) along with the excitation cross section by 
H° impacts computed by Barinovs et al. (2005) 
and electron collision strengths of Wilson &: Bell 
(2002) to compute the C+* energy loss rates ic per 
H atom. These rates are given in the last row of 
Table m They do not change in direct proportion 
to Xe because about half of the excitations come 
from collisions with H atoms, which are more nu- 
merous than the electrons. In addition, the overall 
cooling rate for the WNM does not increase in di- 
rect proportion to ic because it accounts for only 
about one-third of the total cooling. 

Most of the remaining cooling comes from the 
fine-structure excitation of O I and the subsequent 
emissions at 44 and 63 /im. The O I cooling is 
insensitive to changes in Xe because the collision 
rate constants for electrons (Bell et al. 1998) are 
much less than those for neutral hydrogen. Us- 
ing an extrapolation of the atomic hydrogen col- 
lisional rate constants given by Abrahamsson et 
al. (2007) above a temperature of 10^ K, the 
spontaneous decay rates of the two excited lev- 
els given by Galavfs et al. (1997), and 0/H = 
5.75 X 10-"*, we find that for n(H") = 0.47 cm'^ 
and T — 7000 K that the energy loss rate per H 
atom io = 1.13 X lO^^egrg s'^. 

The results for £c shown in the table are 
significantly lower than the average value of 
2 X IQ-^^erg s~^H atom"^ found for low- velocity 



clouds by Lehner et al. (2004) , who measured the 
column densities of C II* from spectra recorded 
by FUSE. [Lehner et al. (2003) obtained simi- 
lar results for sight lines toward WD stars within 
or just outside the Local Bubble.] Lehner et al. 
(2004) compared their results with measurements 
of Ha emission in the same directions, and they 
concluded that about half of the C II* that they 
detected came from fully ionized gas (but with 
large variations from one sight line to the next). 
This inference was supported by the fact that their 
values of £c were slightly lower for sight lines with 
large values of A'^(H I). In principle, one can gain 
an insight on the relative importance of H II re- 
gions by comparing the observed abundances of 
C+* to those of N+*, since the former can come 
from both neutral and ionized regions, while the 
latter arises only from ionized regions. Gry et 
al. (1992) compared these two species in a com- 
prehensive study of both the absorption lines ob- 
served by the Copernicus satellite (Rogerson et 
al. 1973) and the 158 and 205/xm emission lines 
observed by the COBE satellite (Wright et al. 
1991). They concluded that H II regions were re- 
sponsible for a major portion of the C"*"* that was 
observed, but this result is clearly dependent on 
the assumed ratio of atomic C to N in the H II 
regions. 

An important advantage of the tc determina- 
tions synthesized from our results for [Ar I/O I] 
is that they apply only to the WNM; hence they 
give more accurate indications of the carbon cool- 
ing rates within this medium without any contam- 
ination from H II regions. They do, of course, rely 
on the value of assumed relative abundances of C 
and H in the gas phase. 

9.5. Heating Rates 

9.5.1. Thermal Time Constants 

The cooling time for a medium with T = 7000 K 
and thermal pressure uT = 3800 cm^^K (Jenk- 
ins & Tripp 2011) is about 4.1 Myr [ (Wolfire et 
al. 2003); see their Eq. (4)]. However, when the 
medium is impulsively heated and ionized by the 
EUV and X-ray illuminations from a SNR, the 
temperature can approach or exceed 10** K, and 
the onset of La cooling creates a dramatic increase 
in the overall cooling rate (Dalgarno & McCray 
1972). When this happens, the thermal relaxation 
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timescale becomes much shorter than the mean 
interval between the bursts of radiation, each of 
which last only about one to a few times 10^ yr. 
Thus, while we can compute a time-averaged heat- 
ing rate for the SNR illuminations that might ex- 
plain our high levels of partial ionization, we have 
no reason to expect that this average should be 
balanced by the cooling rates that would apply 
for the medium at T = 7000 K. 

9.5.2. Secondary Electrons 

The same energetic electrons that are respon- 
sible for the secondary ionizations Tg and Tg' can 
also heat the gas through collisions with other elec- 
trons. As with the calculations of the efficiency 
of secondary ionizations described in Appendix [XI 
we use the analytic approximations of Ricotti et al. 
(2002) for the heating efficiencies, based on the nu- 
merical results for various conditions obtained by 
Shull & Van Steenberg (1985). Aside from replac- 
ing ionization efficiencies with heating efficiencies 
and multiplying by the energies of the respective 
secondary electrons, the calculations here are vir- 
tually the same as for the ionization rates. 

Table [5] shows the outcomes for the evaluations 
of secondary electron heating rates. It is no sur- 
prise that there is a substantial increase in the 
heating rates when we advance from the ioniza- 
tion created by known sources to either of the two 
hypothetical enhanced ionization examples that 
could explain the electron fractions indicated by 
the observations of [Ar I/O I]. 

Our high result (1.2 x 10~^^erg s~^II atom~^) 
for the heating in the regime where there is in- 
creased penetration of X-rays down to 60 eV 
creates a serious problem for this model, since 
this steady-state rate is considerably larger than 
the corresponding cooling rate £c + — 1.98 x 
10~^® erg s~^II atom~^ given in Section ing 

9.5.3. Cosmic Ray Heating of Electrons 

In addition to ionizing the gas and creating sec- 
ondary electrons that can heat the gas, cosmic 
rays can also interact with free electrons in the 
medium and heat them. The heating rate per unit 
volume is approximately equal to ^CcR,(Il'')n(e), 
where A w 4.6 x lO'^'^erg (Draine 2011, p. 338). 
With the electron densities listed in Table U and 
CcR = 1-25 x lO-i^s-^H atom"\ we find that for 



the three different electron densities in this table 
the heating rates are 1.0 x 10~^^ (for no SNR con- 
tribution), 1.9 X 10"^'' (with SNR contributions), 
and 3.0 x 10~^^erg cm^'^s^^ (for the lower energy 
penetration example). After dividing these num- 
bers by values of n(H°) given in Table lU we find 
that these rates are about 0.4 to 1.2 times the rate 
from secondary electrons liberated by the cosmic 
rays, and 0.05 to 0.5 times the respective heating 
rates from the secondary electrons generated by 
ionizations from the X-ray backgrounds. 

9.6. Dust Grains: Photoelectric Heating 
and Recombination Cooling 

When dust grains are illuminated by starlight, 
they emit photoelectrons, which can heat the 
medium (Watson 1972; Draine 1978; Pottasch et 
al. 1979; Bakes & Tielens 1994; Weingartner 
& Draine 2001b). This heating is partly offset 
by coUisional cooling via grain-ion recombination 
(Draine & Sutin 1987). The efficiencies of these 
mechanisms are regulated by the charge of the 
grains, which in turn depends on n{e), T and G 
(the density of starlight). TablelHshows our calcu- 
lations of grain heating and cooling for the three 
cases (no SNR, with SNR, and lower energy X-ray 
penetration) using the fitting formulae and coeffi- 
cients given by Weingartner & Draine (2001b); see 
note a of the table for details. 

It is generally regarded that the photoelectric 
heating from grains is greater than that from cos- 
mic rays by about one order of magnitude (Draine 
2011, p. 339) [or considerably more than this if 
the assumed value for Ccr(H°) is lower than that 
adopted here, e.g., Wolfire et al. (1995, 2003)]. 
Any enhancement in Xe tends to increase the heat- 
ing rate from grains, since the grains will be more 
negatively charged. However, this increase is not 
as strong as the effect of a greater electron fraction 
on the direct cosmic ray heating, so the disparity 
between the two rates is decreased. 

10. SUMMARY AND CONCLUSIONS 

We have developed a means for deriving the 
representative rates of ionization, and thus the re- 
sulting electron densities, along sight lines that 
penetrate the WNM and that extend out to several 
hundred pc from us, well beyond the edge of the 
Local Bubble. Our method makes use of the fact 
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Table 5 

Secondary Electron Heating Rates' 





Rate 


Source*^ 


(10"27gj.g g-i jj atom"^) 


X-ray background (no SNR) 


4.1 


X-ray background (plus SNR) 


50 


X-ray background (down to 60 eV) 


110 


Main-sequence stars 


1.5^= 


Active binaries 


0.063^= 


WDs 


0.84^= 


Cosmic rays 


4.9^= 


He"*" recomb.'^ (background with no SNR) . 


0.059 


He"*" recomb.'^ (background plus SNR) .... 


1.2 


He"*" recomb.'^ (background down to 60 eV) 


8.8 



^Heating by secondary electrons that are liberated by the ionizations of 
H'^, He'^ and He^, expressed in terms of an energy dissipation rate per 
neutral H atom. Direct interactions of cosmic rays with free electrons 
produce an additional heating which is discussed in Section 19.5.31 

''The conditions for the top and bottom three rows correspond to those 
for the last three columns in Table S) 

'^To first order, the heating rates from internal sources and cosmic rays 
should not change as the strength of the background radiation increases 
above the basic rate from known sources. In reality, they increase by mod- 
est amounts ('--^30%) when the values of Xe (which drive the heating effi- 
ciency) increase. 

'^From secondary electrons that are produced by r//g+(He") and 
r^g+(H'^). Other ionizations arising from helium recombinations do not 
produce electrons with sufficient energy to cause any appreciable heating. 
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Table 6 

Dust Grain Heating and Cooling Rates' 



Without SNR With SNR Greater Low 







Contrib.'^ 




Contrib. 


C 


Enerj 


ry Penetration'^ 


xlO^ 


H 


C 


A 


H 


C 


A 


H 


C 


A 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0.0 


9.9 


5.1 


4.8 


16.1 


7.5 


8.5 


22.0 


10.0 


12.0 


2.0 


15.1 


7.3 


7.8 


23.2 


11.0 


12.1 


30.7 


14.9 


15.8 


4.0 


24.4 


10.6 


13.9 


36.4 


16.0 


20.4 


46.9 


21.8 


25.1 


6.0 


31.6 


13.8 


17.8 


47.4 


20.9 


26.5 


61.3 


28.5 


32.9 



'^Expressed in units of 10~^''erg H atom~^. H = heating, C = coohng, 
and A = H — C (i.e., net heating rate). The values were computed for the 
average interstellar radiation field (ISRF), as defined by Eq. 31 and Table 1 
of Weingartner & Draine (2001b) (with a radiation intensity G = 1.13). The 
adopted value of Ry = Av/E{B — V) — 3.1. The values for the heating and 
cooling rates per H atom tabulated here were computed from Eq. 44 (for grain 
heating) and Eq. 45 (for grain cooling) of Weingartner & Draine (2001b) (but 
without the n(H) factor) and the coefficients listed in their Tables 2 and 3. 

''The abundance of carbon, relative to hydrogen, in the grains. A high value 
of be implies a population of dust grains that is rich in polycyclic aromatic 
hydrocarbons (PAHs). 

'^Thc same conditions as for Columns (2)-(4) in Tabled 
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that when a mostly neutral medium is exposed 
to the ambient EUV and soft X-ray ionizing ra- 
diation, the argon atoms are far more susceptible 
to being ionized than hydrogen atoms. Thus, by 
comparing the abundances of Ar I to those of H I, 
we gain an understanding of the strength of the 
photoionization and secondary processes related 
to it. For a number of practical reasons, we find 
that it is desirable to use O I as a proxy for H I. 
The partial ionization of oxygen is strongly cou- 
pled to that of hydrogen through a rapid charge 
exchange process. Using this strategy, we compare 
the column densities of Ar I and O I derived from 
absorption lines seen in the FUSE spectra of 44 
hot subdwarf stars. 

Since the neutral forms of argon and oxygen 
are virtually absent in regions that are fully ion- 
ized (either the prominent H II regions around hot 
stars or the much lower density but more perva- 
sive WIM), our probes sample only regions that 
have appreciable concentrations of H I. This sets 
our measurements apart from conventional deter- 
minations of average electron densities (e.g., pul- 
sar dispersion measures. Ha intensities, C II fine- 
structure excitation, etc.), which are strongly in- 
fluenced by contributions from the fully ionized 
regions. 

We find that, on average, the abundance of neu- 
tral argon, relative to that of neutral oxygen, is 
[Ar I/O I] = -0.427 ± 0.11 dex below what we 
would expect if both species had no partial ion- 
ization and the solar abundance ratio is a proper 
standard of comparison. We interpret this de- 
ficiency in terms of the greater susceptibility of 
argon to photoionization. After accounting for a 
broad range of processes that can modify the frac- 
tional ionizations, we conclude that with known 
sources of ionizing radiation, both external and 
internal, the only straightforward way to recon- 
cile the large deficiency of Ar I is to propose that 
the neutral medium has a characteristic density 
of only 0.09 cm^^, which is far below a gener- 
ally accepted value of about 0.5 cm~'^. At this 
latter, higher volume density of hydrogen, known 
sources of ionization should produce an electron 
density n{e) = 0.018 cm~^ and create a result 
[Ar I/O I] = —0.14, which is greater than all but a 
small fraction of the measurements and well above 
the overall average value. 

In order to achieve a result for [Ar I/O I] that is 



consistent with both our observations and a den- 
sity n(H I) = 0.5 cm~'^, we must propose noncon- 
ventional sources of additional ionization. We dis- 
cuss two such possibilities; in reality we might be 
witnessing the outcome of some combination of 
both of them working together. 

The first explanation is that the shielding of 
the external, low energy X-ray fixix is less effec- 
tive than expected: if we lower the cutoff energy 
from a nominally expected 90 eV to only 60 eV, 
we can gain sufficient additional ionizing photons 
from the external background to create a level of 
ionization that is consistent with an upper bound 
for the value implied by our observed [Ar I/O I]. 
To explain this lower shielding, one might envi- 
sion a nearly sponge-like topology for the WNM, 
where the internal holes and channels are filled 
with hot. X-ray emitting gas. Indeed, de Avillez 
et al. (2012) have performed hydrodynamical sim- 
ulations of the effects from random SN explosions 
and stellar winds in the plane of the Galaxy, and 
the outcome of their model exhibits a complex, 
turbulent entanglement of hot, warm and cold gas 
complexes. Their simulation, which includes non- 
equilibriimi ionization calculations, yields an av- 
erage electron density n(e) = 0.04 cm~'^, a value 
that is not far removed from our determination 
for the WNM. This near match may be coinciden- 
tal, however. In the simulation, voids containing 
hot, collisionally ionized gas with very low n(e), 
in combination with dense regions that are pho- 
toionized, may produce an outcome that is close 
to the value for the WNM. 

A serious shortcoming of the interpretation that 
we are viewing a steady-state maintenance of a 
high level of ionization in the WNM is that the cal- 
culated heating rate of 1.2x 10~^^ erg s~^ H atom" 
from secondary electrons is unacceptably high - 
much higher than the cooling rate from atomic 
fine-structure excitations or recombinations of ions 
onto dust grains. 

Our second explanation is a time-dependent so- 
lution that can sidestep the problem of overheat- 
ing in the cqiiilibrium case discussed above. Over 
a period of one to a few times 10^ yr, an SNR can 
create a burst of ionizing radiation in the X-ray re- 
gion that can briefly elevate the ionization of the 
ISM to levels well above normal out to distances 
of order a few hundred pc away. Initially, the re- 
combination timescale is quite short, but then it 
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advances to about 1 Myr when the level of ioniza- 
tion approaches our observed average state. For 
heating by secondary electrons, the temperature 
should exhibit a similar rapid adjustment that is 
followed by a much slower trend with time. The 
radiation burst will indeed create a very high heat- 
ing rate, but as the gas temperature approaches or 
exceeds 10^ K, the abrupt onset of the very strong 
La cooling will dump most of the heat over a time 
scale that is considerably shorter than the cooling 
time 4 Myr) for the gas in its usual state. A 
good exposition of this quick thermal recovery to 
T < 10'* K and a comparatively more gentle re- 
laxation in the heat loss and recombinations can 
be seen in plots of T and Xe vs. time in Fig. 2 of 
Gerola et al. (1974). 

We propose that the ionization imprint of an 
SN explosion on the WNM can last well beyond 
the time that the remnant is recognizable in X- 
rays. We have good evidence that a number of 
SNe have exploded in our general vicinity over the 
past 1 — 10 Myr. As a consequence, if one considers 
that the activity in our neighborhood has recently 
been higher than normal for our Galaxy, the WNM 
out to several hundred pc from the Sun may be 
ionized to a level that is somewhat higher than the 
low-density neutral gas at similar galactocentric 
distances elsewhere in the Galactic plane. 
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A. IONIZATION BY SECONDARY ELECTRONS AND Ts 



A.l. Electrons from H and He 

In order to determine rs(H°), [and, later, rs(He'') and rs(Ar'')], we start by considering the effects from 
the secondary electrons that are produced by the primary photoionizations of both H and He. To do so, we 
evaluate various forms of the quantity , Ee,Xe), which represent the number of additional ionizations 
created by the collisions from the photoejected electrons with energy E — IP(H°, Hc+, Hc++) and any 
cascade of additional electrons that may follow. {E is the original photon energy, and IP represents the 
ionization potentials 13.6, 24.6 and 54.4 eV for H and the two stages of He, respectively.) As implied by the 
notation for 0, we must know not only the energy E^_ of the primary electron but also the ambient electron 
fraction Xe = n(e)/[n(H° + n(H+)] since this quantitjO influences the relative efficiency of collisions that can 
create new ionizations, as opposed to collisional heating of free electrons or the excitation of atomic states 
(ShuU 1979). 

For our calculations of (/((H'^, E'e, Xe), we make use of the empirical fits by Ricotti et al. (2002) to the 
secondary ionization efficiencies calculated by ShuU & Van Steenberg (1985)o A complete accounting for 
the secondary ionization of H that includes electrons from the primary ionizations of both H and He is given 
by the relation 



r.(H°) = r,,Ho(HO) + r,3,o(HO) + r,3,+ (H°) 

a(HO,£;)0(H°,i;-IP(HO),:Ee) 



cr(He", £;)0(H°, E - IP(He"), Xe) 
cr(Hc+, £;)0(H", £; - IP(He+), Xe)] F{E) dE 



?i(He°) 
n(HO) 

n(He+) 
n(H") 



(Al) 



Values of cr(He°) and CT(He+) are from Marr & West (1976) and Spitzer (1978, pp. 105-106), respectively. 

Free electrons created by the primary and secondary ionizations of He can exceed 10% of those from H 
and thus should not be neglected. For this reason, we must also consider the helium counterpart to Eq. lAll 
which is given by 

r,,HO (He°) + r,3eO (He°) + r,.He+ (He°) 



r.(He") 



10 

?t(He°) 
n(HO) 
7i(He+) 



i(H") 



cr(H", i;)</>(He°, E - IP(H"), Xe) 
cr(He", E)(j){Re°, E - IP(He°), Xe) 
cr(He+, £;)0(He°, E - IP(He+), Xe)] F{E) dE 



(A2) 



Again, we use the empirical fits of Ricotti et al. (2002) to obtain values of ^(He°, i?e, a^e)- 

In order to compute r3(Ar°), we can estimate from the data compiled by Lennon et al. (1988) and Bell 
et al. (1983) that the cross section for collisional ionization of Ar by low energy electrons arising from H and 



^^ShuU & Van Steenberg (1985) defined their efficiencies in terms of a parameter x = n(H+)/[n(H'') + 7i(H+)] and then assumed 
that the ionization fraction of He is equal to that of H. This assumption is not valid for the conditions considered in the current 
work. Hence, we substitute our term Xe for their x, which allows for a proper accounting of a more generous contribution of 
free electrons from the ionization of He. 

^^For readers who might compare these two references, we provide some guidance about differences in meaning for the represen- 
tations of the efficiencies for ionizations by secondary electrons that each employed: ShuU & Van Steenberg used </> to denote 
an energy fraction that goes into secondary ionizations, whereas Ricotti et al. used $ to mean a number fraction of ionizations 
that take place for each primary ionization. Our function in Egs. lAl] and lA2l is equivalent to <& defined by Ricotti et al. 
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He ionizations is about 3.85 times that for H. The effects of secondary ionizations of Ar° are smaU compared 
to the primary ones, so any moderate deviation from our approximation rs(Ar'^) = 3.85rs(H'^) should be 
inconsequential. 

A. 2. Electrons from Inner Shell Ionizations of Heavy Elements 

In spite of the fact that heavy elements have abundances at least three orders of magnitude below those 
of hydrogen, X-ray ionizations of their inner shell electrons create a perceptible increase in opacity of the 
neutral ISM at energies above 300 eV, and they strongly dominate the absorption of X-rays above about 
1 keV (Morrison & McCammon 1983; Wilms et al. 2000). These ionizations create energetic electrons, which, 
like those from the primary ionizations of H and He, can create additional collisional ionizations Tg' of the 
species that we are considering. 

Electrons from the ionization of the inner shells of atoms arise from two fundamental processes. First, an 
electron is ejected by the primary ionization, and it has an energy equal to that of the ionizing photon minus 
the ionization potential of the inner shell electron. To calculate the ionization cross section ij{Xi_j^E) for 
level J of an element Xi as a function of photon energy E, we used the coefficients and fitting formula given by 
Verner & Yakovlev (1995). The primary ionization is followed by the ejection of one or more Auger electrons, 
which can also collisionally ionize H, He, and Ar. The energies Ei^j of the Auger electrons, or the sum of 
energies of two or more electrons, have been calculated by Kaastra & Mewe (1993). In our assessments of 
the strengths of the resulting collisional ionizations of H, He, and Ar, we use equations analogous to Egs lAll 
and Ell 

r,, (H") = {A{X^) [ . ^) [pi^^^E - IP(X,j), xe) + (biR",E,^j,Xe)\ F{E) dE^ (A3) 

and 

r,,(He°) = 10^{a(X,)^ j a{X,^j,E)\cp{Y{e\E~W{X,^,),Xe) + m(^\E,^3^x,)\F{E)dE'^ , (A4) 

i 3 

where A{Xi) is the assumed abundance of element X relative to that of H and the (f> functions are identical 
to those described in Section lA.il As we did for the secondary ionizations of Ar, we assumed a rate equal to 
that of H multiplied by a factor of 3.85. For the Auger electrons, we have no information on the distribution 
function of electrons with energy. Thus, we fall back on the simple assumption that the overall effect is the 
same as that for a single Auger electron with an energy equal to the average sum of energies of multiple 
electrons (if more than one). In doing so, we overlook the shortcoming that when two or more electrons have 
almost equal energies, the efficiency for secondary ionizations should be slightly lower than that for a single 
electron with the same energy. 

For the element abundances Ax, we assume that only the gas-phase (depleted) abundances are relevant. 
While it is true that high energy X-rays can penetrate dust grains and ionize the inner shells of electrons 
for atoms inside the grains, the escape probabilities for these electrons are low. For instance, Weingartner 
et al. (2006) calculated that dust grains with radii a — O.l /im have ultimate electron yields less than 30% 
per photoionization for E < 1 keV. In various models for the grain size distributions, about half of the mass 
of material is in grains with a either above or below about 0.1 /xm (Draine 2011, p. 281). For the gas-phase 
depletions, we adopt values given by Jenkins (2009) for low density gas, i.e., strengths characterized by his 
parameter = 0. The elements considered were, in decreasing order of importance, O, C, Ne, N, S, Si, and 
Mg. These elements accounted for more than 99% of the free heavy atoms in the ISM. 

Finally, we follow the strategy of Adamkovics et al. (2011) in neglecting ionizations caused by fluorescence 
photons from heavy elements that are exposed to X-rays. They showed that this is an effect that is small 
compared to those from primary and Auger electrons. 



43 



B. IONIZATION BY PHOTONS FROM THE RECOMBINATION OF HELIUM IONS 



When a free electron with an energy of order kT combines with a He ion and creates a lower ionization 
stage of He (with an ionization potential IP) in either a level with an energy i?o above the ground electronic 
state or in the ground state itself (£'0 = 0), a photon having an energy E = W — Eq + kT is generated. 
A transition between bound levels also emits a photon with an energy E equal to the energy difference of 
the levels. Since the mean free path of a photon that can ionize H or He is short compared to the size of 
the neutral gas region, it is safe to assume that any recombination radiation photon capable of creating an 
additional ionization will do so and not escape. (We neglect the absorption of photons by dust grains.) 

In the following subsections, we consider the various ways that the photons from He recombinations can 
be produced and how they add to the overall ionization rates. In most instances, they can ionize more than 
one kind of atom or ion. Thus, as with the treatment of primary ionizations by internal sources discussed in 
Section VTM we need to consider the relative probability y{X,E) defined earlier by Eq. [^Hthat any photon 
of energy E will ionize a given species X. This sharing function will appear in the detailed treatments that 
follow. 

While the only the rates Tfje+O^'^) and rjLfgo(H") (and their counterparts for Ar) are of direct interest in 
the development of the relative ionizations of H and Ar, the influence of helium recombination photons in 
ionizing He° and He+ will affect the ionization balance equations given in Sections 16.41 and 16.51 Except for 
the relatively rare recombinations to the He"*" ground Is level that can then photoionize H*^ (see Section lB.3l 
below), ionizations of H and He arising from He recombinations do not create free electrons with enough 
energy to contribute to Tg for either H or He. 

B.l. rj,g+(He+) from the Reaction He++ + e -)• He+ + 7 

The only recombination channel that creates a photon with an energy that is capable of reionizing He+ is 
that which goes directly from the free electron continuum to the Is ground state of He"*". This recombination 
results in a photoionization rate 



with T4 — T/W^K (Draine 2011, p. 138; all future equations of this form are from the same source). Were 
it not for the fact that some fraction of the emitted photons ultimately ionize He" and H°, as indicated by 
the fact that the y term in Eq. IBll is less than 1, the effect of r^g+(He^) would be simply to convert the 
familiar Case A recombination rate (to all electronic levels) to that of a Case B rate (to levels with n > 2) 
(Baker & Menzel 1938). 

B.2. rjje+ (He°) from the Reaction He++ + e ^ He+ + 7 

As with the ionization of He+, some of the photons that arise from recombinations directly to the Is state 
will ionize He", giving a rate equal to that expressed in Eas. lBll and lB21 except that the term ?/(He", 54.4eV+ 
kT) replaces y(Re'^ , 54.4eV + fcr). These photons are supplemented by direct recombinations to the 2p level, 
which then decays to the ground state by emitting a 40.8 eV photon. Finally, recombinations to levels above 
n = 2 eventually emit photons with an average energy of about 50 eV. (This average energy can degrade 
toward 40.8 eV if 7T.(He~'')/n(H") is large enough to create resonant absorptions of the He+ Lyman series 
photons before they can ionize H".) These three processes yield an ionization rate 



rffe+(Hc+) = l"^''^ ai,(He+, T);/(He+, 54.4eV + kT) s^^ , 




(Bl) 



where an approximate fit for the recombination coefficient to the Is state is given by the formula 



ai,(He+,r) - 3.16 x io-i3(T4/4)-"-5«-"-°i^'"(^-/4) cm^s'^ , 



(B2) 



rffe+(He°) 



n(He++)n(e) 
n(He") 
+ [aB(He+,T)- 



{ais(He+, r)y(He°, 54.4 eV + kT) + a2p(He+, T)2/(He°, 40.8 eV) 
- a2p(He+,T) - aeff.2s(He+,r)]2/(He°,50eV)} , 



(B3) 
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where the direct rate to just the 2p state is given by 

a2p(He+,T) = 1.07 x io-i3(T4/4)-"-68i-"-06iin(T4/4) ^^3g-i ^ 

recombinations to all levels except the ground Is state have a combined rate constant 

aB(He+, T) = 5.18 X 10-13(T4/4)-0-833-0.0351n(T4/4) ^^3g-l ^ 

and all recombinations that pass through the 2s state have a combined rate 

acff,2s(He+,r) = 1.68 X 10-13(r4/4)-°-™5-0.0081 ln(T4/4) ^^3g-l^ (Bg) 

[The last equation is a Draine-style fit to the three numbers in Tabic 14.2 of Draine (2011).] 

B.3. rjje+(H°) from the Reaction He++ + e ^ He+ + (1 or 2)7 

There are multiple processes that generate photons capable of ionizing H'^. First, there are direct recom- 
binations to the Is state. Direct recombinations to the 2p state, followed by decays to the Is level, first 
generate a photon with an energy E = 13.6 eV + kT that can ionize hydrogen, and then produce a 40.8 eV 
photon that can ionize either H" or He°. Recombinations to the metastable 2s state again start with a 
creation of a hydrogen-ionizing photon, and later the decay of this level creates 2 photons, which have a 
combined average hydrogen ionization efficiency of 1.42 photons (Osterbrock 1989, p. 29). As discussed 
earlier in connection with the derivation of THe+ (He"), photons with « 50 eV are emitted by decays from 
levels with n > 2. Finally, there are cascades from upper levels to the 2s state that generate only 1.42 
photons without the added effect of the initial E ~ 13.6 eV -I- kT photon. For the processes that have just 
been discussed, we have 

rHe+ (H") = ^^^^|J^|^{ai.(He+, r)y(H°, 54.4eV + kT) + a2p(He+, r)[l + y(H", 40.8 eV)] 

+ 2.42a2«(He+,T) + [aB(He+,T) -a2p(He+,T) -ae//,2«(He+,T)]y(H°,50eV) 

+ 1.42[a,;y,2s(He+,T)-a2.(He+,r)]} s-i , (B7) 

where 

a2.(He+,T) = 4.68 x io-i4(y^/4)-o.537-o.oi9i„(T4/4) ^gg) 

B.4. rjje+(Ar+) from the Reaction He++ + e He+ + 7 

The ionization potential of Ar+ is 27.62 eV, which is only about 3 eV higher than that of He'^. There are 
no photons for cascades of levels within recently formed He+ that have energies that fall between the two 
ionization potentials, hence Ea. lB3l can be used for the ionization of Ar+, after substituting every occurrence 
of "HeO" with "Ar+." 

B.5. rjje+(Ar°) from the Reaction He++ + e He+ + (1 or 2)7 

The processes that ionize Ar are very similar to those expressed for H in Section IB. 31 except that direct 
decays from the continuum to the 2p or 2s level do not have enough energy to ionize Ar'^. The absence of 
these two sources of photons for ionizing Ar" eliminates several terms that we had in Eq. IB7[ which now 
leads to 

r^,+ (Ar°) = ^(He++Me) {^^^(He+, r)y(Ar°, 54.4 eV + kT) + a2p(He+, r)y(Ar°, 40.8 eV) 
n(Ar ) 

-I- [ae (He+ , T) - a2p(He+ ,T)-a,ff^ 2. (He+ , T)] y (Ar*' , 50 eV) 

+ 1.29a,/^,2.(Hc+,r)[n(Ar")/n(H°)][(a(Ar"))^(^)/(a(H°))^.(^)]}s-^ , (B9) 
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For Ar°, the factor 1.42 for the effective number of H ionizing photons arising from the 2-photon decay of 
the 2s level must be reduced because the ionization potential of Ar" (15.76 eV) is slightly higher than that of 
H° (13.60 eV). After integrating the photon energy distribution f{E) given by Spitzer & Greenstein (1951) 
from the ionization potential of Ar" to 40.8 eV, we find that the photon yield is reduced to a value equal to 
1.29. Since these photons are spread over a broad range of energy, we must calculate a flux-weighted ratio 
of the photoionization cross sections of Ar and H that appears in Eq. IB9[ given by 



<.40.8cV / <.40.8cV 

(a(Ar°))^.(^)/(a(H°))/(s) = / f{E)a{kv^)dE / f{E)a{B°)dE, (BIO) 

JlP(ArO) / JIP(HO) 



which has a numerical value of 13.8. 

B.6. ThbO (He°) from the Reaction He+ + e ^ He° + 7 

As with the case for He+ reionizing itself discussed in Section [B.l) we have for He" 

rHeo(HcO) = "^^''^y ai.(He",T)y(He", 24.6 eV + kT) s'^ , (BU) 
n(Hc ) 

where an approximate fit for the recombination coefficient to the Is state is given by the formula 

Qi,(He°, T) = 1.54 X 10-13^-0.486 ^^^3^-1 ^3^2) 

(Benjamin et al. 1999). 

B.7. r/jeo(H°) from the Reaction He+ + e He° + 7 

There are two channels that contribute to rjifeo(H"). The first is direct decays to the Is state. The second 
is all possible decays to levels with n > 2, of which 96% of them generate photons which can ionize H° in 
the low density limit (Osterbrock 1989, p. 26). Hence, we find that 

r^e«(H°) = ^^^^^^^-^^[ai,(He°,r)y(H°,24.6eV + fcr)+0.96aB(He°,T)] s^^ (B13) 

where 

aB(He",r) = 2.72 x lO-^Tf"-^^^ cm^s"^ (B14) 

(Benjamin et al. 1999). 

B.8. r//eo(Ar°) from the Reaction He+ + e ^ He° + 7 

The ionization of Ar*^ from recombinations to He'^ is similar to that for H" given in Eq. IB13I 

r«eo(Ar°) = »(He+)n(e) ^^^^(^^o^ T)y{Kr\ 24.6 eV + kT) 
71 (Ar ) 

+ 0.94aB(He°,T)n(Ar")CT(Ar", 20eV)/[n(H")o-(H", 20eV)]} s"^ (B15) 

Drake et al. (1969) have computed the energy distribution of the two-photon decay from the singlet 2s level 
of He'^. From this distribution, one finds that the average yield of H- ionizing photons from each decay is 
0.56. With the slightly higher ionization potential of Ar", this yield drops to 0.30 photons, which lowers the 
coefficient in front of aB(He°) from 0.96 to 0.94. The average energy of all these decays is 20 eV, which must 
be used in the calculation of the ratio of the Ar primary photoionization cross section to that of H. 
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C. COSMIC RAY IONIZATION 



We consider that a cosmic ray ionization rate for molecular hydrogen has a representative rate Ccr(H2) = 
3.5 X 10~^^s~^ (Indriolo & McCall 2012). (There are large variations found from one region to the next 
however.) This rate includes both primary cosmic ray ionizations and the additional ionizations from the 
resulting secondary electrons, and it should be twice as large as that for atomic hydrogen (Nagy & Vegh 
1992). The measurements of Indriolo & McCall (2012) applied to cool, diffuse cloud regiona^i where some 
molecules are present. Since these regions have a^e ~ 3 x 10^^ (Draine 2011, p. 185), we must adjust (qr 
to reflect the fact that the efhciency of secondary ionizations is reduced when larger values of Xe are present 
(see Appendix E)) . According to Spitzer (1978, p. 144) the average energy of such secondary electrons is 
35 eV. Neglecting the fact that there is some spread in the actual energies, we derive an ionization rate for 
primary ionizations alone for H in atomic form, 

Cp,cr(H") - 2[l + (/.(H0,35cV,3x 10-4)] ' ^^^^ 

(See Section [A. II for an explanation of the term (/)(H°, i?e, Xe).) Within the context of the far more diffuse 
regions that are studied here, which have much higher electron fractions Xg, we must allow for a reduction 
in the ionization efficiency for the secondary electrons to obtain a corrected total rate for both primary and 
secondary ionizations, 

Ccr(H°) = [1 + 0(HO, 35 eV, a;e)]Cp,CR(H") . (C2) 

As with their susceptibility to ionization by photons, argon atoms are about 10 times more easily ionized 
by protons than hydrogen (for proton energies E ^ 1 MeV) (Hooper et al. 1962; Kingston 1965). In the 
description of secondary ionizations in Section lA.ll it was noted that at low electron energies the cross 
section for coUisional ionization for Ar is about 3.85 times that of hydrogen. Accordingly, we find that 

CcR(Ar°) = [10 + 3.850(H", 35 eV, Xe)]Cp,CR(H°) . (C3) 

Finally, as one of the ingredients for calculating P'(Ar) (Section 16. 5|) and the He recombination photons 
(Section [B]) involves the ionization equilibrium of He fEgs. [TUl -ll5 p . we must know in principle the cosmic- 
ray ionization rate of He, 

CcR(He") = 1.36Cp,cr(H") (C4) 

Again, we used the results of Hooper et al. (1962) and Kingston (1965) to find that Cp,CR(He°)/Cp,CR(H°) = 
1.36. Most of the secondary electrons are not energetic enough to ionize He, so Eg. IC4I does not have a 
factor included. In practice, CcR(He°) ^ r(He°). 



^^These regions are well shielded from low energy X-rays, and hence virtually all of the ionization is caused by cosmic rays. Of 
the 50 sight lines covered in the survey by Indriolo & McCall (2012), only one (e Per) had N(H I -|- 2H2) < 10^^ cm~^. Under 
these circumstances the X-ray ionization rate should be less than 10~^* i.e., far less than the rate from cosmic rays. 
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